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12. ANÁLOGOS APLICÁVEIS A PROBLEMAS 
DE CONDUÇÃO DE CALOR 


Os problemas de condução de calor podem, 
em geral, reduzir-se à resolução de equações do 
tipo 


T 
1) div (kgrad T)=cs: : R 
dt 
sendo 
t o tempo 
T a temperatura 
c o calor específico 
c a densidade 
k tensor condutividade térmica, que se reduz 


a um escalar no caso de meios isótropos. 


No caso de sistemas isótropos uniformes com 
k independente de T, a equação 1) reduz-se a 


d 
2) lap T= Ea: 
x dt 
sendo 
3) dam cê , 


Consideremos um circuito eléctrico constituído 
por um reticulado de resistências, R, ligando nós, 
e capacidades ligando esses nós a um nó comum, 
O, de referência das tensões. 
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Dentro da aproximação de diferenças finitas, 
as tensões dos nós, u, em relação ao nó O, satis- 
fazem à equação 


4) div (7 gradu)=C — ; 


ot 


sendo 7 o tensor condutividade (em determinada 
«escala» geométrica) do reticulado e C a capa- 
cidade por unidade de volume dos nós (na mes- 
ma escala geométrica). 


Estabelecendo a correspondência, 


t o» É 
| 
g mem js O 
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E Ce. 
ú Es» E 


corrente <. > fluxo de calor, 


tem-se, imediatamente, um análogo eléctrico do 
sistema térmico traduzido por 1). 

Suponhamos, por exemplo, o caso de uma 
tubagem metálica isolada, com simetria cilíndrica, 
conforme representado na fig. 12-1, inicialmente 
a temperatura uniforme, e que se pretende deter- 
minar o regime transitório subsequente a um 
aquecimento brusco, por exemplo por introdução 
de vapor cuja temperatura seja mantida pratica- 
mente constante. 
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Fig. 12-x — Representação esquemática de uma tubagem metálica isolada. 
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O problema tem simetria cilíndrica e axial, 
bastando considerar em troço de comprimento 
41 e ângulo 12, conforme representado na 
fig. 12-2. 


Fig. 12-2 — Troço da tubagem da fig. 12-1 considerado 
para análise de um problema de condução de calor 
com simetria cilíndrica e axial, 


Decomponhamos este troço, conforme indicado 
na fig. 12-3, em várias camadas, definidas por 
superfícies S,. Seja R, a resistência térmica (radial) 
entre a superfície i a superfície i+1 e seja C, a 
semi-soma das capacidades térmicas dos troços 
definidos pelas superfícies S, ,eS,eS,eS,.. 

Ter-se-á o análogo indicado na fig. 12-4, sendo 
R,e €C, no análogo, capacidades e resistências 
proporcionais a R e €,, no sistema térmico. 

A simulação do regime transitório correspon- 
derá a aplicar uma tensão escalão u, (análogo da 
aplicação brusca de uma temperatura T, à super- 
fície S,). A evolução no tempo das tensões nos 
vários nós é o análogo da evolução no tempo das 
temperaturas nas várias camadas S (fig. 12-5). 

Suponhamos agora o caso de uma fonte de 
calor, por exemplo destinada a aquecimento da 


Fig. 12-3 — Decomposição do troço de tubagem da fig. 

12-2 em camadas definidas pelas supefícies S,, para 

estudo de um problema de transmissão de calor com 
simetria cilíndrica e axial. 


tubagem, correspondente à introdução de calor 
em secções do tubagem igualmente espaçadas, con- 
forme representado esquematicamente na fig. 12-6. 
Trata-se de um problema redutível a duas dimen- 


——— mm 


t 


Fig. 12-5 — Representação esquemática da evolução no 
tempo das tensões u,, análogas das temperaturas das 
superfícies S,, correspondente a uma variação brusca 
da temperatura da superfície S, (tendo por análogou, ), 
a partir dz uma situação de equilíbrio térmico. por 
exemplo por admissão brusca de vapor na tubagem. 


Fig. 12-4 — Análogo eléctrico do troço de tubagem da fig. 12-3, para efeitos de condução de calor. A resis- 

tência r representa as condições de transmissão de calor ao exterior (supostas equivalentes a condução se r 

for constante); quando tal aproximação não seja válida, r será substituído por um elemento não linear, com 
corrente, por exemplo função não linear de u. 
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Fig. 12-6 — Representação esquemática de uma tubagem metálica isolada com anéis de aquecimento 
circulares uniformemente espaçados. 


sões considerando um sector de tubagem de ângulo 
] “ A . 

x e comprimento = é sendo 1 a distância entre 

anéis de aquecimento consecutivos. 


Fazendo uma decomposição desse sector, con- 
forme esquematizado na fig. 12-7, tem-se um 


Fig. 12-7 — Sector da tubagem da fig. 12-6 considerado 
para efeitos de condução de calor. Atendendo às 
condições de simetria da tubagem, para causas de 
aquecimento com igual simetria basta considerar um 


| 
troço de ângulo 4x e comprimento ——, sendo 1 a 
2 


distância entre os anéis de aquecimento. 
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análogo do tipo do representado na fig. 12-8, 
correspondendo as resistência às resistências tér- 
micas entre facetas paralelas do sector segundo 


. , R 1 
o raio e segundo o eixo, e as capacidades a Z da 


soma das capacidades caloríficas dos troços do 
sector vizinhos da aresta do sector (definida por 
um valor de i e um valor de j) correspondente a 
cada nó do análogo (definido por um valor de i 
e um valor de j). Haverá também que considerar 
a transmissão de calor ao exterior, de tempera- 
tura mantida constante, o que, a título ilustrativo, 
se supõe corresponder à condução representada 
pelas resistências r. 


A introdução brusca nos anéis de aquecimento 
de uma potência térmica constante corresponde 
a injecção no nó do modelo correspondente ao 
anel (definido por j=1 e i=-3) de uma corrente 
escalão (amáloga de metade da potência térmica 
no anel, no ângulo Àz, por o análogo correspon- 
der a metade da zona de influência de cada 
anel) conforme representado na fig. 12-8. A evo- 
lução no tempo das tensões nos vários nós 
é o análogo da evolução mo tempo das tempera- 
turas dos pontos correspondentes a esses nós 
(fig. 12-9). 

O tipo de análogo anteriormente descrito usa 
resistências e capacidades. Outra solução, é con- 
siderar, em vez das capacidades, resistências 
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Fig. 12-8 — Análogo elécrrico do troço de tubagem 
representado na fig. 12-7. O nó O representa o nó de 
referência de tensões; para comodidade de represen- 
tação, não é representada a ligação entre os vários 
terminais ligados a este nó, que são todos referenciados 
por O, As resistências r representam as condições de 
transmissão de calor ao exterior, supostas equivalentes 
a condução se r for constante ; quando tal aproximação 
não seja válida, r será substituída por um elemento 
não linear, com corrente, por exemplo, função não 
linear da tensão aos seus terminais. À referenciação 
dos nós e do circuito, segundo os índices i e j. corres- 


ponde às intersecções das superfícies Se s,, com a 
notação da fig. 12-7, e a tensão do nó ij, em relação ao 
nó de referência, O, é o análogo da temperatura nessa 
intersecção. A corrente i, injectada no nó 31, é o aná- 
logo de metade da potência térmica transmitida à 
tubagem por um sector 54 de um anel de aquecimento. 
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Fig. 12-g — Representação esquemática, para uma po- 
tência de aquecimento variando segundo uma função 
escalão, a partir de uma situação de equilíbrio, da 
evolução no tempo da corrente i (correspondente à 
potência térmica de aquecimento), e das tensões em 
vários pontos (correspondentes a temperaturas), para 
o análogo da fig. 12-8, correspondente a uma estilização, 
representada na fig. 12-7, do comportamento térmico 
da tubagem indicada na fig. 12-6. 


ligadas a uma fonte de tensão, como represen- 
tado esquematicamente na fig. 12-10 para o caso 
de um sistema a uma dimensão, fazendo o cál- 
culo e ajuste das correntes injectadas nos nós 
para uma sucessão de instantes discretos, 


) t=(0;1;...;n— 1;n;n do DO vu bri, 


Fig. 12-10 — Análogo eléctrico para um problema de transmissão de calor a uma dimensão, para resolução 
considerando sucessivos instantes de tempo e ajuste em cada uma das correntes i. As resistências r, os 
potenciómetros -, e a fonte de tensão (contínua ou alternada de amplitude constante) destinam-se a obter 
as correntes i que simulam as correntes nos condensadores segundo um análogo do tipo da de fig. 12-4. 
Para referência das tensões dos nós, poderá adoptar-se, por exemplo o nó O, obtido de forma indicada na figura, 
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ajustando a corrente i; injectada num nó i no 
instante n At pela relação 


6) f=C(uj—u, ), 


ou outra semelhante, sendo C, a «capacidade» 
que, no modelo com resistências e capacidades, 
seria ligada no nó i. 

A vantagem deste tipo de analogia é usar 
apenas resistências, sendo o análogo construtiva- 
mente mais simples. Em contrapartida, há neces- 
sidade de ajuste das correntes e, por vezes, de 
considerar um número de instantes n At muito 
elevado, o que torna o método moroso, se não 
for parcialmente automatizado. Há calculadores 
baseados neste método em que é efectuado 
automáticamente o ajuste das correntes k (1). 


(1) Um exemplo de calculador deste tipo é tratado por 
]J. 1. Missen e D.W. Stops no artigo The solution of 
heat transfer problems by automated Liebmann net- 
works, Proceedings of the International Association 
for Analog Computation, Vol. XI, n.º 1, de Janeiro 
de 1969. 
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C. D. U. 624.341.5/94 


ELEMENTOS SOBRE A PRODUÇÃO E O CONSUMO DE ENERGIA 
NA REDE ELÉCTRICA NACIONAL **) 


I — Breve nota mensal 


As afluências ao conjunto do sistema produtor atingi- 
ram cerca de 1186 GWh, valor que tem uma probabi- 
lidade de ser excedida “da ordem dos 25", e a que 
corresponde um coeficiente de produtibilidade hidroe- 
léctrica de 1,33. 

As afluências registadas desde o início do ano 
hidrológico (trimestre Novembro-Janeiro) tém uma 
probabilidade de serem excedidas da ordem dos 27'º/,. 

Para efeito da continuação de um regime de afluén- 
cias favorável, durante o mês de Janeiro o armazena- 
mento no conjunto das albufeiras registou um acréscimo 
de 419 GWh, vindo a atingir no fim do més cerca de 
1592 GWh, valor que corresponde a 75º/, do armaze- 
namento máximo possível. 


LI — Elementos gerais (GWh) 


a) Moss 


o Variação 
1967 | 1968] “lu 
Produção hidráulica (Ph)... «| 371,7, 590,5 + 59 


Produção térmica (Pt)..... “| 981) 17,5|— 82 
Produção total (Hr). | 469,8 | 608,0 / + vy 


Energia recebida de empresas 


não pertencentes ao RNU (Er) 1,1 5,0 +409 
Exportações (Ex) ,..ec,.c ce. | 14 2,6 |— 86 
Importações (1).... ... Sea ca 48,5 4,0|— 92 
Saldo importador (Sp)....... +46,9 | +IAI— 97 
Cons, em bomb, hidroelé, Ch) (1) 0,0 29,9 | 


Produção para con- 
sumos perman. (Pep)..... 901,1 
Produção para con- 
sumos não perman, (Pcnp) -».» 
DOtãl vestes sos xa 
Coeficiente de hidraulicidade 
—————e— DDD — 
NOTA: 
(1) O aumento percentual da produção para consumos permanen- 


tes, tendo em conta a incidência dos domingos e dias especiais 
t 6.3. 


(1) 
5823 + 6,2 


| 16,7 | 528/4216 
917,8 | 615,0 + 19 
0,84] 1,83] — 


II — Diagramas de carga dos dias característicos 


| 4.º teira: 


Produção hidráulica 
Produção térmica 
Produção total 


(E n) MWB /12001 
(Pro)MWh | 3882 1350 
(Pr) MWb [15883 |2º311 


Trocas com f Export. (Ex) MWh oz no 
Espanha ( Import.  ( 1)MWh | 227% 70 
Cons. em bomb, hidroel, (Cp) MWh U 1390 


Prod. p.* cons. perm.  (Pep) MWh 17 506 
Prod.p.*cons.não perm.(Penp)MWh 600 1856 


TOTAL Pr-+(-Ex) MWh | 18 106 


E | 


Potência máx. MW. 1400 1147 
Potência min. MW 414 509 


Pr + (Ex) 


:E Utiliz. da ponta horas 18,1 18,2 
E E Factor de carga 0,75 0,76 
PR Potência máx. MW To | 1061 
ES p Potência min. MW 389 426 
a SP | Utiliz. da ponta horas 18,0 17,9 
Factor de carga 0,75 0,75 
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19687 | 1968 


JANEIRO 


2.3 - DisGRAMAS DE CARGA DOS DIAS CAJACTERÍSTICOS (Pro) Es-C,) 
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MR PRN armazenada nas principais albufeiras 


No fim do mês 
Albufeiras: E, ES 

GWh 0 13) 
Alto Rabagão . .....c| 3689 51,9 
PaEEÓOIA . é ds =» 65 & =| 128 26,2 
Venda Nova . « «. «cc ««| 1098 85,8 
Salamonde +. +. +. +... E a 14,8 03,6 
Caniçada É cd E CS IO, or jm 14,3 43,2 
Vilar . cce es re] 59,1 58,7 
Cabril so ss cs cc vaj MA 81,8 
Castelo do Bode +... ... 96,9 59,2 
Guilhofrei . «. . cc cv co 7,5 91,6 
Lagoa Comprida « . +... 34,4 (2) 92,9 
Hanta LUMA - ve we. é 45,5 73,9 
Pracana » cc cvs c ici 11,3 87,5 
PO xa Es uk é ga 8,70) 65,4 
e, com A, Rabagão . . .| 1173,7 55,8 
Potal | sem A, Rabagão. . «| “804,8 70,0 


NOTAS 
(1) Coeficiente de enchimento em energia das albufeiras, 


(2) Loclui 1,9 GWh armazenados em Vale do Rossim no início do 
més e 2,4 GWh no fim do més, 


(3) Inclui 1,9 (GWh armazenados no açude do Poio no início do mês 
e 2,7 GWh no fim do mês. 


(*) Elementos extraídos das estatísticas mensais do Repar- 
tidor Nacional de Cargas (R. N. C.). As produções e 
vos consumos das empresas do R.N.C. representam 
cerca de 94, dos totais do Pais. 
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ANÁLISE ELÁSTICA DE VALES E TALUDES 


RESUMO 


Assimila-se a crusta terrestre a um semi-espaço elás- 
tico onde se instalam estados planos de deformação. 
Formula-se assim um problema de elasticidade plana 
relativo a um semi-plano indefinido, 

Considera-se a resolução deste problema por dois 
métodos numéricos, o método dos elementos finitos e um 
metodo integral, ecomparam-seos resultados obtidos para 
um caso concreto, o da determinação das tensões num 
estrato argiloso que ocorre no talude das Portas do Sol 
em Santarém. 


1 — DEFINIÇÃO DO PROBLEMA 


Consagra-se o presente trabalho à determina- 
ção de tensões originadas pela abertura de vales 
na crustra terrestre. 

Não se discute a legitimidade do modelo mate- 
mático que servirá para o cálculo dessas tensões: 
não se fará mais que apresentar e comparar dois 
métodos numéricos aproximados para a resolução 
de um problema de elasticidade plana que se 
admite, sem apresentar justificação, como repre- 
sentação razoável da realidade. 

Assimila-se a crustra terrestre a um semi- 
-espaço elástico em que se instalam estados planos 
de deformação. Reduz-se assim o problema, que 
em rigor é tridimensional, a um problema bidi- 
mensional relativo a um semi-plano indefinido 
homogéneo e isótropo. 

Considere-se um sistema cartesiano de eixos 
ortogonais e suponha-se que o semi-plano corres- 
ponde a x, > O (fig. 1). 


&2 


Fig. 1 


O campo de tensões instalado deve equilibrar 
as forças exteriores. Não contando com as 
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SYNOPSIS 


By identifying the earth crust with an elastic half- 
-space subjected to a plane state of strain, the elastic 
analysis ofvalleys and slopes can be converted into a plane 
elasticity problem defined on a half-plane., 

The author discusses the solution of this problem by 
two numerical methods, the fnite element method and an 
integral method, and compares the results yielded by both 
when applied to an actual case: the andlysis ofstresses in 
a clay layer which occurs at the Portas do Sol slope in 
Santarém. 


reacções aplicadas no infinito, essas forças são 


unicamente as devidas à acção da gravidade, 
Resumem-se assim a forças de massa de densi- 


dade 
X a 7i, (2) 


sendo o peso específico do material e i, o versor 
associado com o eixo Ox,. 

Supondo que a fronteira do semi-plano se 
encontra livre de tensões, e atendendo a que, em 
virtude da simetria elástica e geométrica, o campo 
de tensões instalado não deve depender de x,, a 
forma analítica mais geral que as componentes 
de tensão podem apresentar, sem que se viole 
nem o equilíbrio nem a compatibilidade, é 


Gp=ax+b (3) 
É ag = — JX, (4) 
p= 0 (5) 


ae b são constantes cujo valor não poderá ser 
determinado sem formular hipóteses de ordem 
geológica. Considerações deste tipo permitem 
admitir como razoável uma tensão 9, propor- 
cional a 7, isto é, tomar b = 0. Fazendo 
a ==— «y vem pois 


= To (6) 


No caso de um líquido, « seria igual à unidade. 
A experiência mostra que o valor de « é em 


regra mais baixo, da ordem do coeficiente de 
Poisson. 

Ao abrir-se um vale de perfil ABCD (fig. 2), 
o campo de tensões inicial é perturbado na res- 
pectiva vizinhança, já que o troço ABCD fica 
livre de tensões aplicadas. 


Fig. 2 


Para determinar o novo campo de tensões, 
pode-se raciocinar da seguinte maneira: 

Considere-se novamente o semi-plano indefi- 
nido e calculem-se as componentes 7, e 7,, rela- 
tivas aos eixos coordenados, das tensões aplicadas 
ao longo de um troço de linha, ABCD, corres- 
pondente ao perfil do vale (fig. 3). Sendo v o 
angulo da normal à linha com o eixo 0x,, tem-se: 


7) = 7,98eN P— 0, COS P (8) 


Fig. 3 


Seja agora o semi-plano com o entalhe corres- 
pondente ao vale (fig. 2). O campo perturbador 
é o que se desenvolve neste domínio se se apli- 
carem no troço ABCD, como forças exteriores, 
tensões de componentes iguais e de sinal contrário 
às definidas por (7) e (8), mantendo-se nulas as 
tensões em pontos situados a distância infinita 
do entalhe (fig. 4). 
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Com efeito, se se sobrepuser o campo elástico 
correspondente a estas forças ao instalado inicial- 
mente no semi-plano, obtém-se um campo elástico 
equilibrado pelas forças de massa definidas pela 
equação (1), ficando livre de tensões o perfil do 
vale e a fronteira do semi-plano para um e outro 
lado do entalhe. 

Por outro lado, as tensões aplicadas no infinito 
não são afectadas pela abertura do vale, já que 
as do campo perturbador tendem a anular-se. 

O problema elástico que este trabalho se propõe 
resolver é o da determinação do campo pertur- 
bador. 

Tal problema não tem resolução formal conhe- 
cida para um perfil de vale inteiramente arbitrário. 

A análise do campo elástico introduzido por 
uma força concentrada na fronteira do semi-plano 
indefinido (problema de Flamant) revelar-se-á 
contudo da maior utilidade. 

Por esta razão se lhe dedica o próximo parágrafo. 


2— O PROBLEMA DE FLAMANT 


Considere-se um sistema de forças paralelas, 


com resultante F, distribuídas numa região do 
semi-plano indefinido a que pertença um ponto P 
da fronteira, e suponha-se que a região se contrai 
progressivamente até se reduzir ao ponto P, 
mantendo-se porém constante a intensidade da 
força. 


No limite tem-se uma força concentrada F 
actuando na fronteira do semi-plano (fig. 5). O 
problema de Flamant consiste na determinação 
do campo elástico correspondente e este limite. 


Fig. 5 


Uma primeira observação a fazer é a de que, 
podem obter-se forças concentradas actuando no 
plano indefinido, como limites de processos dife- 
rentes do que foi descrito. É por exemplo o caso 
de um entalhe semicircular de raio R, praticado 
na fronteira do semi-plano, nos bordos do qual 
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se distribuem uniformemente forças radiais de 
resultante F (fig. 6). Supondo que R tende para 


zero, mantendo-se F constante e uniforme a 
distribuição, encontra-se no limite outro campo 
que só a uma distância infinita do entalhe se 
confunde com o de Flamant, tanto no que se 
refere a tensões como a deslocamentos. 


Fig. 6 


Uma segunda observação é que o problema 
não fica perfeitamente definido desde que não 
se introduza a condição de que as tensões se 
anulam a uma distância infinita do ponto de 
aplicação da força. 

Na realidade a fronteira total do domínio não 
é so a recta x,=0, que normalmente recebe o 
nome de fronteira, mas também a fronteira no 


infinito. A força Ti equilibrada por forças apli- 
cadas nessa fronteira, e é importante reconhecer 
que as componentes da densidade de distribuição 
destas forças, isto é, as tensões aplicadas no 
infinito, podem ser nulas, mesmo que a resultante 
o não seja. 

Não interessa discutir em pormenor o problema 
de Flamant. Basta saber que as expressões das 
tensões e dos deslocamentos são as seguintes: 


2F cos(9— 0) 


é + r (9) 
9 =9,4=0 (10) 
u, = — EA (9-0) log see 
x E. 
—- TÊ E (0) sen (8-8) (11) 
E 
2u F “mu 2 É r 
u,= 22 E sen(-1)+ 2 log E sen (1-9)— 
E ” ô 
; PA. 


1— e 
es . E H (0-5) cos (0-0) + E F sen (1-0) 
(12) 
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A fig. 7 define as grandezas que intervêm 
nestas fórmulas. 


Fig. 7 


O campo de deslocamentos definido pelas 
expressões (11) e (12) é distorsional, e não pode- 
ria portanto aparecer como solução no caso de o 
problema dizer respeito ao plano indefinido total, 
e não só a um sector do plano. 

Vê-se que só as tensões radiais são diferentes 
de zero. Esta propriedade é característica dos 
campos gerados por forças actuando no vértice 
de cunhas (o serni-plano é uma cunha cujo ângulo 
é de 180º), e uma distribuição de tensões que a 
apresente diz-se distribuição radial simples. 

Outra propriedade interessante consiste em o 
campo elástico não depender da orientação da 
fronteira do plano mas somente da da força, o 
que se conclui por as grandezas dependerem de 0 
através de 9-0... 

O ponto P é um ponto singular do campo: 
tensões e deslocamentos apresentam nesse ponto 
valores infinitos. Fora desse ponto, o andamento 
das tensões radiais é perfeitamente razoável: decres- 
cem ao longo dos raios vectores, até que se anulam 
no infinito. 

Já a variação dos deslocamentos é paradoxal e 
merece um comentário mais atento. É que a fun- 
ção log r, de que dependem linearmente as com- 
ponentes, se torna infinita não só para r=? como 
para 1=9, 

Esta particularidade poderá parecer estranha, 
especialmente se se atender a que no problema 
de Boussinesq, isto é, no do semi-espaço (tridi- 
mensional) solicitado por uma força concentrada, 
tanto tensões como deslocamentos se anulam no 
infinito. 

Não convém porém esquecer que os problemas 
de deformações planas são na realidade proble- 
mas tridimensionais em que as forças actuam 
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desde x,= — co a x,=+oc(fig. 8). A força F 
é somente a que aplica na porção do semi-espaço 
compreendida entre dois planos perpendiculares 
a x, e distanciados de uma unidade. A força total 
aplicada é portanto infinita, e isso explica que se 
obtenham no infinito condições diferentes das do 
problema de Boussinesq. É evidente que tais 
condições não se apresentam na Natureza. 


Fig. 8 


Uma consequência prática é a de que, enquanto 
no problema de Boussinesq está indicado eliminar 
os deslocamentos de corpo rígido atribuindo o 
valor zero aos deslocamentos no infinito, no pro- 
blema de Flamant há que fixar arbitráriamente a 
origem dos deslocamentos. 

As fórmulas (11) e (12) foram assim escritas 
na hipótese de os deslocamentos radiais se anu- 
larem em um ponto situado sobre a recta de acção 
da força a uma distância r, do seu ponto de 
aplicação, a qual é perfeitamente arbitrária. A 
rotação do semi-plano foi porém escolhida de 
maneira que os pontos colocados sobre a recta 
de acção da força permanecem sobre a mesma 
recta após a deformação. 

É interessante notar que, para 0—9 =—+, 
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Quer dizer: os pontos situados sobre a recta 
que passa por P e é perpendicular à direcção da 
força têm deslocamento radial constante para a 
direita de P e para a esquerda de P, desde r=0 
a r=o0, Em P, o deslocamento radial sofre uma 
discontinuidade: mantém-se o valor absoluto, 
mas altera-se o sinal, de forma que todos os 
pontos tendem a aproximar-se de P. 

No problema que se pretende resolver, as 
forças que actuam no entalhe têm uma resultante 


R diferente de zero, igual ao peso do material 
eliminado. A uma distância do entalhe suficiente- 
mente grande, em comparação com as dimensões 
deste, o campo perturbador confunde-se com o 


campo gerado por uma força concentrada R 
actuando no semi-plano indefinido (fig. 9). 


Fig. 9 


O comportamento dos deslocamentos e das 
tensões a essa distância é assim o mesmo que no 
caso do problema de Flamant. 

Isto permite reduzir o problema, que foi 
definido para um domínio infinito, a um de domí- 
nio finito. 

Basta supor o domínio semi-infinito cortado 
ao longo de uma linha EFGHI (fig. 9). A acção 
da parte do semi-plano exterior a essa linha pode 
ser substituída, de acordo com a hipótese funda- 
mental admitida na Teoria da Elasticidade, por 
forças actuando ao longo de EFGHI. 

Essas forças poderiam ser determinadas com 
toda a facilidade desde que se conhecesse o 
campo de tensões instalado no domínio com 
entalhe. Na prática, não se conhece, pois são 
justamente essas tensões que se deseja determinar. 
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Sabe-se porém que, se a fronteira EFGHI esti- 
ver suficientemente afastada do entalhe, a distri- 
buição de tensões se confunde com a devida à 


força R suposta concentrada na fronteira de semi- 
-plano. Isto permite estabelecer condições de 
fronteira ao longo de EFGHI, definindo-se assim o 
problema para um domínio finito. 

Em vez de tensões, podem fixar-se desloca- 
mentos. Esses deslocamentos poderão ser também 
os determinados no problema de Flamant. 


3 — MÉTODOS NUMÉRICOS NA ELASTICI- 
DADE PLANA 


Chegou a altura de se indicar quais os métodos 
numéricos de que se dispõe para a resolução de 
problemas de elasticidade plana. Supõem-se nulas 
as forças da massa. 

Trata-se portanto de determinar soluções numé- 
ricas aproximadas para o seguinte conjunto de 
equações de campo: 


equações de equilíbrio : 


Cit ot X = 
491 too TX =0 


relações deformações-deslocamentos : 


“o ="/2 (u,9 Hu5,) 


relações tensões-deformações (deformações planas, 
meios isótropos) : 


E 
“uv q+)) (1-2) À dio 
Prado in 
2 qro a—2 (UT meu 
pre + 
Lo q+y! 


Estas equações de campo devem verificar duas 
condições por ponto da fronteira, as quais podem 
envolver tensões ou deslocamentos. Tem-se assim 
o primeiro ou o segundo problemas fundamentais 
da elasticidade plana. 

Dispõe-se de três espécies de métodos aproxi- 
mados para resolver este problema: o método 
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das diferenças finitas, o método dos elementos 
finitos e os métodos integrais. 

Nos dois primeiros, a análise numérica começa 
antes da integração das equações de campo. Nos 
últimos, a análise numérica parte de integrais 
gerais das equações de campo que se particula- 
rizam para as condições de fronteira do problema, 
as quais só aproximadamente são satisfeitas. 

Em qualquer caso, o fim a atingir é discretizar 
o problema contínuo, isto é, fazê-lo depender de 
um número finito de parâmetros cuja determina- 
ção se faz resolvendo um sistema de equações 
lineares. 

Não se considera neste trabalho o método das 
diferenças finitas. Trata-se de um método que 
parece substituído com vantagens pelo método 
dos elementos finitos, mais fácil de automatizar 
e susceptível de interpretações físicas mais 
directas. 

O parágrafo seguinte dedicar-se-á à exposição 
do método dos elementos finitos aplicado à elas- 
ticidade plana. 

Focar-se-ão depois os métodos integrais e mos- 
trar-se-á como podem ser especialmente adapta- 
dos para o tratamento do problema dos vales e 
taludes. 


4 — MÉTODO DOS ELEMENTOS FINITOS 


Considere-se o domínio D, bidimensional, para 
o qual se define um determinado problema de 
elasticidade plana. Suponhamos este domínio divi- 
dido em um número finito de subdomínios e, 
dentro de cada subdomínio, famílias de campos 
dependentes de um número finito de parâmetros. 

Um elemento finito é cada um destes sub- 
domínios juntamente com a família dos campos 
nele admissíveis. Qualquer campo desta família 
é uma combinação linear com coeficientes q; de 
um número finito de modos fundamentais de 
deformação. Cada campo da família corresponde 
a atribuir valores particulares aos parâmetros q. 

Para parâmetros, são em regra escolhidas as 
componentes de deslocamento em pontos situados 
nas fronteiras dos elementos, os quais se dizem 
nós ou pontos nodais. Os deslocamentos nodais 
são os deslocamentos generalizados. 

O tipo de um dado elemento depende da sua 
forma genérica, da especificação dos seus pontos 
nodais e dos campos admissíveis, que se definem 
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analiticamente exprimindo o campo de desloca- 


Ed 


mentos, u, em termos dos parâmetros q; . Sendo 
u o vector das componentes de deslocamento e 
q o vector dos parâmetros, vem: 
=f(x,, x) 4 (17) 
Os elementos normalmente utilizados são trian- 
gulares ou quadrangulares. Destes últimos só os 
rectangulares serão considerados neste trabalho. 
O tipo de elemento triangular mais frequente- 
mente empregado caracteriza-se por se tomarem 
os vértices para nós e por se supor linear a 
variação dos deslocamentos : 


| ui = Ap + Ay X + AX 


| U = Ay + Ay X + AX o 

A matriz f da equação (17) pode ser determi- 
nada a partir da equação (18). Com efeito, desig- 
nando por 1, Il e III os vêrtices do triângulo, e 
por x, ; “A coordenadas 
(fig. 10), vem: 


as respectivas 


mW ti 
HI (x, ,x2) 


(x) x!) 
HE (od 


Fig. 10 


= Aw+ Am x, + Ai X 


u! 
| As q 
u, = Aw+ Ay x + Ap x, 
Em virtude de ser 
q=|u, 
ul 
; (20) 
I 


vem, mais compactamente, 


q=MA (21) 
em que: 
M-[1 84 0 & 071 
Die q ama 
LA mM a 6 O 
(22) 
a. 
HI HI 
1 x” xy 0 O 0 | 
| 
HI HI 
106 60 1% 
e 
= ( , 
A Ai 
A, 
A 
12 
half (23) 
Am 
Ag 
( As 


A matriz M é uma matriz quadrada não-singular 
(desde que os vértices do triângulo não sejam 
colineares). Poder-se-á portanto explicitar A em 
função de q: 

A=M'q (24) 

Por outro lado, a equação (12) pode tomar a 
seguinte forma matricial 


= XA (25) 
sendo 
Jim OQ 
k = | (26) 
| O O O 1 x, X, 
Introduza-se (24) em (25). Resulta: 
u=XM-!q (27) 


Comparando (27) com (17) conclui-se final- 
mente ; 
TAN (28) 


Para os elementos rectangulares é frequente 
utilizar-se a distribuição de deslocamentos 
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uU,=Aj + AX + AjX + AjsX,X 
tg SE Ay + AgX, + An + Aos XX 


Os diferentes elementos são compatibilizados 
por condições reduzidas de compatibilidade. 
Estas impõem a continuidade das componentes 
de deslocamento em nós coincidentes de ele- 
mentos adjacentes. 

Supondo que os lados dos rectângulos são 
paralelos aos eixos coordenados, conclui-se das 
expressões (29) que permanecem rectos após a 
deformação. O mesmo acontece com os lados 
dos triângulos. 

Nestas condições, as distribuições de desloca- 
mentos (18) e (29) têm a vantagem de permitir 
que os lados coincidentes de elementos adjacentes 
se mantenham colados após a deformação, isto é, 
que as condições reduzidas de compatibilidade 
sejam suficientes para garantir a continuidade 
dos deslocamentos em todo o domínio, e, por- 
tanto, a compatibilidade. 

A energia de deformação E” de um elemento 
pode ser determinada considerando a equação (17). 

Com efeito, se as relações deformações-desloca- 
mentos (15) tomarem a forma simbólica 


e=Du (30) 


em que D é um operador diferencial e 


“n 
e = | *w (31) 
28,9 


e as relações tensões-deformações (16) a forma 
(32) 


em que 


(33) 
sendo H uma matriz simétrica, vem, para expres- 
são de densidade da energia de deformação: 
T T T 
W=!/s e=!,0e He=!- (Du) H(Du)= 
T F 
= qo (Df) H(Df) q 


(34) 
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Logo 
T 
Eº = [w dS=4hp q* Kºqf (35) 
ge 


em que 


E 
Kº = fr  H(D f) ds (36) 
Je 


A matriz Kº diz-se matriz de rigidez do ele- 
mento e. 

Recorde-se agora o teorema da energia poten- 
cial total minima que estabelece que a solução 
exacta do problema elástico é, de todas as solu- 
ções compatíveis, a que torna mínima a energia 
potencial total, P. 

Compreende-se que a solução aproximada 
gerada pelos métodos dos elementos finitos seja 
a que minimiza P, não na classe de todas as 
soluções compatíveis, mas somente na dos campos 
admissíveis. 

A energia potencial total é a soma da energia 
de deformação com a energia potencial das forças 
exteriores: 


P=E+P, (37) 


Tem-se 
T 


E=2 ES =1/s 2 q Kº q" (38) 


e, em virtude de serem nulas as forças de massa, 


n T 
P=— | p u dc (38) 


“1 


em que p representa as forças exteriores actuando 
na porção c, da fronteira do corpo (onde foram 
fixadas tensões). O traço superior indica tratar-se 
de uma grandeza fixada. 

Tem-se pois: 


E Tr B: ; 
Pp md pp ude=—3 [9 f de qó (39) 
e e 
cê ce 


Neste somatório serão diferentes de zero sô- 
mente os termos correspondentes aos elementos 
que tiverem pontos na porção c, da fronteira. 

Diz-se vector das forças generalizadas actuando 
no elemento & o vector 


4 à 
q = J p fdc (40) 
ce 
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À energia potencial das forças exteriores pode 
pois tomar a forma : 


T 
P, = 0º qó (41) 


A expressão de P é assim: 
T 2 
p=3[ th q Kº q — Qº | (42) 
e 


Em virtude das condições reduzidas de compa- 
tibilidade, os vectores qº não são todos indepen- 
dentes: os deslocamentos dos vértices coincidentes 
de um grupo de elementos adjacentes têm de ser 
iguais. 

Diz-se nó do sistema de elementos um ponto 
que seja nó de um ou mais elementos. Os deslo- 
camentos nodais do sistema são componentes de 
deslocamento desses pontos. Seja q o respectivo 
vector, 

As condições reduzidas de compatibilidade são 
expressas por relações lineares entre os vectores 


q e o vector q: 
q =Tq (43) 


Substituindo em (42), vem: 


T T 
= q Kg —Q q eo 
em que 
T q 
K=5Tº KCTº (45) 
T A 
Q=5Tº 0º (46) 


Como se disse, a aproximação introduzida pelo 
método dos elementos finitos consiste em mini- 
mizar a energia potencial total, não na classe de 
todas as soluções compatíveis, mas na dos campos 
admissíveis, definidos pelas matrizes f* e pelas 
relações (43). 

Dado que os elementos do vector q são todos 
independentes, as condições de estacionaridade 
podem ser determinadas diferenciando sucessiva- 
mente a equação (44) em ordem a todos os ele- 
mentos de q e igualando a zero. Resulta a equação 
matricial ; 
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Kkq= 0 (47) 


No caso de os deslocamentos nodais fixados 
serem em número suficiente para eliminar os des- 
locamentos rígidos, tornando o vector q determi- 
nado, podem-se extrair da equação (47) os 
deslocamentos nodais desconhecidos. 

Os campos de deslocamentos no interior de 
cada elemento podem determinar-se a partir do 
vector Q, associando as equações (43) e (17). 
Vem: 


u=fTº q no elemento e. (48) 
Conhecido o campo de deslocamentos, podem 
determinar-se as tensões por: 


s=HDu (49) 


Tal é na sua essência o método dos elementos 
finitos. 

É importante assinalar que as constantes elás- 
ticas podem variar de elemento para elemento, 
pelo que podem ser tratados, com toda a facili- 
dade, meios heterogéneos. 


5 — MÉTODOS INTEGRAIS 


Podem incluir-se nos métodos integrais os 
métodos baseados na resolução de equações inte- 
grais e na representação conforme do domínio 
num círculo de raio unitário. 

O método que vai ser apresentado é de natu- 
reza diferente e tem sido descrito em vários 
trabalhos [1]. 

Os métodos integrais, conforme se afirmou, 
caracterizam-se por conduzirem a soluções que 
satisfazem exactamente as equações de campo e 
baseiam-se no princípio da sobreposição, tomando 
como solução do problema uma combinação linear 
de N soluções elementares, linearmente indepen- 
dentes, das equações de campo, multiplicadas por 
coeficientes que são determinados obrigando a 
solução a verificar N condições em diferentes 
pontos da fronteira. 

É-se assim naturalmente levado a um sistema 
de equações lineares. Contudo, enquanto no 
método dos elementos finitos as incógnitas (os 
deslocamentos nodais) são pontos localizados 
tanto na fronteira como no interior do domínio, 
nos métodos integrais a ordem do sistema é 
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geralmente muito mais baixa, já que as incógnitas 
se referem sômente aos pontos da fronteira. 


Esta circunstância, que permite uma utilização 
mais eficiente de capacidade de memorização do 
computador, é compensada pelo facto de a 
matriz do sistema conter, no método dos ele- 
mentos finitos, muitos elementos nulos, o que não 
acontece nos métodos integrais. 


Se nos métodos integrais se considerar um 
número infinito de soluções elementares, o sis- 
tema de equações transforma-se numa equação 
integral. Tal acontece nos métodos de Miche- 
-Massonnet e de Sherman-Lauricella. 


No método que vai ser descrito isto não acon- 
tece. 

Na realidade, embora uma equação integral 
seja muito útil na teoria, sê-lo-á muito menos na 
análise numérica, não só porque as passagens ao 
limite restringem frequentemente a generalidade 
dos métodos (a isso se deve por exemplo a difi- 
culdade no tratamento de problemas mistos), mas 
também porque o rigor que introduzem é só 
aparente, já que raramente é possível determinar 
soluções formais destas equações. 


A primeira característica do método que se 
apresenta é considerar um número finito de solu- 
ções elementares. A sua primeira vantagem, a 
capacidade para tratar problemas mistos. 


Por outro lado, os métodos de Massonnet e 
Sherman não podem ser aplicados na resolução 
de problemas definidos para domínios múltiplica- 
mente conexos, pois a solução aproximada que 
constróem resulta da sobreposição de campos 
distorsionais, que são na realidade os correspon- 
dentes à distribuição radial simples. 


As soluções elementares no método em questão 
correspondem a campos uniformes de desloca- 
mentos, isto é, não são distorsionais, podendo o 
domínio ser multiplamente conexo. 


O nosso método pode assim resolver os três 
problemas fundamentais para domínios simples 
ou múltiplamente conexos. 


Seja um domínio D, simples ou múltiplamente 
conexo, e as correspondentes condições de fron- 
teira. Deseja-se determinar o campo elástico no 
interior de D. 

Suponha-se o domínio imerso no plano inde- 
finido ou, mais precisamente, suponha-se que os 
contornos fechados que formam a fronteira são 
desenhados no plano indefinido. 
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Estes contornos limitam uma área R+ do plano 
com a mesma forma que o domínio. 

Seja R a área complementar e F a fronteira 
comum às duas áreas. 

Considerem-se agora contornos auxiliares à 
volta de cada um dos contornos da fronteira 
efectiva. Os contornos auxiliares são exteriores a 
R” e seu conjunto forma a fronteira auxiliar T. 

Considerem-se ainda forças auxiliares actuando 
no plano indefinido, distribuídas sobre T, de 
forma que as condições de fronteira do domínio 
D sejam reproduzidas em F. Se tal distribuição 
existir, o campo que ela gera na área RT é o 
campo procurado pois: 


1 — verifica as equações de campo em todos os 
pontos (pois que as forças auxiliares se 
distribuem em R”, de forma que R” fica 
livre de forças exteriores e singularidades). 

2 — as condições de fronteira são verificadas 
na fronteira. 


Assim, se puder ser determinada uma distri- 
buição de forças em T tal que as condições de 
fronteira sejam reproduzidas em F, o campo em 
D fica conhecido, visto que sabemos sempre deter- 
minar, em qualquer ponto do plano indefinido, 
os efeitos de forças concentradas noutros pontos. 

Existirá porém sempre uma distribuição em T 
de forças fictícias (forças auxiliares) que possam 
reproduzir em F as condições de frónteira de D? 
Se esta distribuição existir, será ela única? É fácil 
mostrar que a condição para a existência de tal 
distribuição é que o campo em D possa ser 
prolongado analiticamente para fora de F, de 
forma a estender-se a todos os pontos de R-. 
Demonstra-se também que, se existir, será única. 

Em geral, a prolongação analítica até grande 
distância da fronteira não é possível, pelo que 
há que empregar contornos auxiliares muito 
próximos dos contornos efectivos. O sistema de 
equações vem então, por outro lado, muito melhor 
condicionado. Em contrapartida a reprodução das 
condições de fronteira é fortemente influenciada 
pelo andamento poligonal da fronteira auxiliar 
e do diagrama das forças auxiliares. 

A discretização pode obter-se desde que se 
considere a distribuição de forças na fronteira T 
como função de um número finito de parâmetros. 
Quando esse número tender para infinito, obter- 
-se-á a verificação das condições de fronteira em 
todos os pontos desta. Na prática, só se pode 
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conseguir uma verificação das condições de 
fronteira em um número finito de pontos. 

A distribuição de forças poderia ser definida 
de muitas maneiras. A mais simples parece 
consistir em tomar tantos pontos de cada contorno 
auxiliar (pontos auxiliares) quantos os pontos 
dados em cada contorno efectivo para os quais 
haja que respeitar as condições de fronteira. 
Seguidamente, não há mais que supor que a 
densidade de distribuição varia linearmente entre 
cada dois pontos e que o contorno auxiliar é 
poligonal com vértices nos mesmos pontos. 


A distribuição depende assim de 2 N parâmetros, 
os N pares de componentes da densidade das 
forças auxiliares nos N pontos auxiliares. 

Cada solução elementar refere-se ao valor 
unitário de cada uma das incógnitas, tomando o 
valor zero todas as restantes. 

Há assim que resolver o seguinte problema 
elementar : determinar as tensões e deslocamentos 
em cada ponto do plano indefinido devidas a 
uma distribuição de forças paralelas, actuando 
em segmentos AA” e AA“ com uma extremidade 
comum. 

A distribuição supõe-se linear em ambos os 
segmentos, anulando-se a densidade em A” e A” 
e tomando em À o valor zero. 

Este problema elementar é facilmente resolvido, 
integrando os deslocamentos e tensões devidos 
a uma força concentrada actuando no plano 
indefinido. 

O programa que foi preparado, de acordo com 
estes princípios, resolve os três problemas funda- 
mentais para meios homogéneos e isótropos 
(poderiam ser anisótropos), domínios simples ou 
multiplamente conexos, finitos ou infinitos de 3.º 
espécie, na ausência de forças de massa. 


O programa pode ser utilizado de duas maneiras. 
Na primeira, são dadas condições de fronteira, e 
o computador calcula e imprime as componentes 
das densidades das forças auxiliares nos pontos 
auxiliares e, a partir dessas componentes, calcula 
tensões e deslocamentos em qualquer ponto do 
domínio. 

Deve notar-se que as condições de fronteira 
podem ser deslocamentos ou forças aplicadas em 
troços contíguos da fronteira, nunca tensões. De 
facto, se se fixam as intensidades das resultantes 
das tensões em troços contíguos, e sendo os troços 
suficientemente pequenos para se fixarem apro- 
ximadamente as posições dessas resultantes, o 
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princípio de Saint-Venant garante que o campo 
obtido seja praticamente coincidente com o campo 
efectivo a alguma distância do contorno, depen- 
dente da distância entre pontos. O mesmo não 
aconteceria se se fixassem tensões. 

Na segunda maneira de se utilizar o programa, 
as densidades das forças auxiliares são directa- 
mente dadas e determinam-se tensões e desloca- 
mentos em qualquer ponto. 

Isto permite, por um lado, não executar de 
uma só vez todos os cálculos: desde que se 
guardem as forças auxiliares que se obtiveram 
da primeira vez, podem determinar-se tensões e 
deslocamentos noutros pontos, em qualquer altura. 

Por outro lado, podem determinar-se efeitos 
de forças que actuam em linhas do plano inde- 
finido, desde que se assimilem tais linhas a 
poligonais, e se considere também como poligonal 
o diagrama da densidade de distribuição das 
forças que actuam nessas linhas. 

Esta última facilidade permite que o programa 
seja por exemplo aplicado na determinação de 
efeitos de forças de massas, etc. 


6—UTILIZAÇÃO DOS MÉTODOS 
NUMÉRICOS NA ANÁLISE DE VALES 
E TALUDES 


A determinação do campo perturbador constitui 
um problema de elasticidade plana de tipo especial. 

O domínio é por um lado simplesmente conexo, 
por outro é infinito. 

Viu-se já que o problema pode ser reduzido 
ao da análise de um domínio finito. Tal é o 
caminho mais óbvio e será esse em geral o caminho 
a seguir se se quiser aplicar o método dos elementos 
finitos. 

No caso do método integral, não se obtiveram 
bons resultados quando se quis transformar o 
domínio num domínio finito utilizando uma 
fronteira artificial. 

A razão do insucesso foi muito simplesmente 
a de que os pontos a considerar na fronteira têm 
de ser suficientemente próximos uns dos outros 
para que as condições de fronteira possam ser 
convenientemente reproduzidas. Dada a extensão 
da fronteira, o número de pontos a empregar 
resultou demasiadamente grande para a capacidade 
do computador. 

Nestas condições houve que pensar numa 
adaptação do método, sacrificando embora a sua 
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generalidade, que não seria necessária nas presen- 
tes condições. 

A primeira observação que se fez foi a de que, 
sendo o domínio simplesmente conexo, não have- 
ria inconveniente em empregar como soluções ele- 
mentares distribuições radiais simples, desde que 
a semi-recta de descontinuidade dos deslocamentos 
destes campos ficasse situada fora do domínio 
(fig. 11). 


SEMI-RECTAS DE DESCONTINUIDADE 


Fig. 11 


Considere-se um semi-plano em cuja fronteira 
actua uma distribuição de forças paralelas a uma 
certa direcção do plano e cuja densidade, p, tem 
uma distribuição poligonal com vértices em um 
certo número de pontos definidos na referida 
fronteira (fig. 12). 


A / 


Fig. 12 


O campo gerado por esta distribuição pode ser 
calculado se se tiver calculado o campo gerado 
por um diagrama triangular com base em três 
pontos arbitrários, A, B e €, da fronteira do 
semi-plano (fig. 13): não há mais que considerar 
uma distribuição unitária com base em cada ponto 
da distribuição e nos dois pontos vizinhos e 
multiplicá-la pela correspondente ordenada do 
diagrama. Com efeito a fronteira do plano, para 
a esquerda do ponto A e para a direita do ponto 
B, ficará livre de tensões. 


AB T 


Fig. 13 
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Considere-se agora o semi-plano com um enta- 
lhe e forças actuando no entalhe (fig. 14) segundo 
um diagrama poligonal. Considerem-se da mesma 
forma diagramas triangulares centrados em vários 
pontos considerados no entalhe e sobreponham-se 
esses diagramas multiplicados pelas correspon- 
dentes ordenadas do diagrama. Agora, a fronteira 
do semi-plano, para a esquerda do ponto C e 
para a direita do ponto D, já não fica livre de 
tensões. 


Fig. 14 


No nosso método consideramos portanto uma 
distribuição de forças auxiliares desde um ponto 
E a um ponto F (fig. 15), escolhida de forma 
que as tensões se anulem ao longo dos troços 
EC e DF e tomem no entalhe os valores fixados. 


Fig. 15 


A distribuição é poligonal, dependendo por- 
tanto de um número finito de parâmetros e a 
condições de anulamento das tensões só é res- 
peitada aproximadamente, impondo o anulamento 
das forças aplicadas a troços sucessivos da fron- 
teira. 

Para a esquerda de E e para a direita de F, as 
tensões aplicadas poderão ser consideradas nulas 
se os troços EC e DF foram muito grandes em 
comparação com as dimensões características do 
entalhe. 

Considere-se com efeito (fig. 16) uma distri- 


Fig. 16 
buição radial simples centrada num ponto P a 
uma distância à da fronteira. Para pontos situa- 


dos na fronteira, a uma distância | de P, as 
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componentes de tensão 7, e 7,, aplicadas na fron- 
teira são: 


q 
tendendo portanto para zero com T Se a força 


for perpendicular à fronteira, este efeito é refor- 
çado pelo facto de q, ser proporcional a sen x. 

Podem portanto considerar-se como nulas as 
tensões aplicadas na fronteira do semi-plano e 
tomar para solução aproximada a gerada pela 
distribuição de forças entre E e F. 

Para adaptar o programa do método integral 
ao problema dos vales e taludes, a modificação 
essencial consistiu assim em empregar como 
campos elementares distribuições radiais simples, 
e não, como no programa geral, campos gerados 
por forças concentradas no plano indefinido. 

Por outro lado a fronteira auxiliar limita-se a 
cobrir a zona EF da fronteira real, sendo portanto 
finita sem ser fechada. 


7 — EXEMPLO DE APLICAÇÃO 


Os princípios que acabaram de ser expostos 
foram aplicados na resolução de um problema 
posto à D.M.A. pelo Serviço de Geotecnia. 
Tratava-se da determinação de tensões num 
estrato argiloso, à profundidade de 73 metros, 
que ocorre na encosta das Portas do Sol em 
Santarém. Tal determinação revestiu-se de parti- 
cular interesse para a concepção e dimensiona- 
mento das obras que aí foram executadas. 


Embora o problema se refira, como foi dito, a 
um talude, preferiu-se considerar um vale simé- 
trico (fig. 17), o que, se por um lado, como é 
intuitivo, não altera sensivelmente o valor das 
tensões no estrato a investigar, por outro permitiu 
introduzir uma simetria que simplificou grande- 


mente a análise numérica. 


Fig. 17 


Tomou-se para o peso específico das terras, y, 
o valor de 2 t/mº, adoptando-se para 2 o valor 0,5. 
As forças exteriores correspondentes ao problema 
perturbador estão esquematizadas na fig. 18. 


ps 200 tim? 
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Fig. 18 


O método dos elementos finitos foi aplicado 
mediante a decomposição da estrutura (em vir- 
tude da simetria, pôde ser considerada somente 
metade da estrutura total) em 102 elementos (57 
rectângulos e 55 triângulos) a que correspondia 
102 nós (figs. 20 e 21). Levou assim a um sistema 
de 204 equações. Sendo 35 a largura da semi-banda 
de elementos significativos, tem-se, atendendo 
à simetria da matriz de rigidez, necessidade de 
memorizar pelo menos 6720 coeficientes. 

O método integral exigiu a consideração de 
63 pontos na fronteira do domínio (fig. 19). Em 
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Fig. 19 — Método integral — Pontos da fronteira 
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Fig. 20 — Decomposição em elementos. Conjunto geral 


virtude de existir um eixo de simetria, a ordem 
do sistema reduziu-se também a 63, o que, 
tratando-se de um sistema não simétrico em que 
todos os coeficientes são significativos, implica 
a memorização de 3969 elementos. 


100 


300 


No que se refere a tempos, o método integral 
exigiu 200 minutos (no computador National- 
-Elliot 803-B) até ao fim da resolução do sistema, 
mais 4 minutos por cada ponto onde se procedeu 
à determinação de tensões e deslocamentos. Com 


Fig. 21 — Decomposição em elementos. Pormenor 
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o método dos elementos finitos, gastaram-se 35 
minutos na formação da matriz de rigidez, 205 
na resolução do sistema (cujas incógnitas são os 
deslocamentos) e 35 na determinação das tensões 
a partir dos deslocamentos : ao todo 275 minutos. 

O método dos elementos finitos fornece os 
deslocamentos em todos os nós e as tensões nos 
centros de todos os elementos. Se, pelo método 


20 há 


À 100 


Metodo dos elementos finitos e 


Metodo integral 


| 


+ 


Tensões G, co longo da recta AB 


integral, se tivessem determinado tensões e des- 
locamentos em 20 pontos, gastar-se-iam 280 mi- 
nutos, o que praticamente seria o mesmo tempo 
que no caso do método dos elementos finitos. 
Uma exploração mais pormenorizada do campo 
levaria a maiores tempos de cálculo. 

As figuras 22 a 27 apresentam os resultados 
obtidos. 


Fig. 22 — Tensões «,; ao longo do estrato AB 


8 — CONCLUSÕES 


Tanto o método dos elementos finitos como o 
método integral foram aplicados na resolução de 
um mesmo problema. 

O método dos elementos finitos implicou a 
memorização de maior número de coeficientes. 
No que se refere a tempos de cálculo, pode-se 
porém considerar que os dois métodos se equi- 
valem, já que em regra o conhecimento das ten- 
sões e deslocamentos em 20 pontos (tinha-se 
verificado que, para este número de pontos, o 
tempo correspondente ao método integral era o 
mesmo que o correspondente ao dos elementos 
finitos) constitui informação suficiente. 

O andamento das tensões fornecidas pelo mé- 
todo integral apresenta-se, por outro lado, mais 
razoável que o das fornecidas pelo dos elementos 
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finitos : neste último caso, com efeito, as tensões 
apresentam flutuações pronunciadas nas regiões 
cobertas por triângulos. Tais flutuações não têm 
no entanto importância de maior, e a conclusão 
a tirar é a de que convém decompor em rectân- 
gulos as zonas onde se pretender informação 
mais precisa sobre as tensões, pois que os ele- 
mentos rectangulares, dado que consentem varia- 
ção linear das tensões, conduzem a melhores 
resultados. 

Na prática adoptar-se-ão rectângulos sempre 
que possível, reservando os triângulos para onde 
não puderem dispensar-se, isto é, para a vizi- 
nhança da fronteira e para a transição entre ele- 
mentos rectangulares de diferentes tamanhos. 

Os elementos de comparação seriam diferentes 
se diferente número de elementos finitos ou de 
pontos da fronteira tivesse sido considerados 
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em cada um dos métodos. O exemplo apreciado 
parece no entanto típico, pelo que será lícito afir- 
mar que o método integral conduz, para os mes- 
mos tempos de cálculo, a diagramas de tensões 
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Metodo dos elementos finitos » 


Metodo integral 
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com melhor aspecto e a melhor utilização da 
capacidade da máquina. 

Ambos os métodos são no entanto satisfatórios 
e não convém esquecer as vantagens do método 
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Fig. 23 — Tensõe c,9 ao longo do estrato AB 
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Metodo dos elementos finitos a 


Metodo integral 
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Tensões C,pio longo da recta AB 


Fig. 24 — Tensões C,, ao longo do estrato AB 


dos elementos finitos: consente heterogeneidades 
do terreno, a consideração de furos ou aberturas 
e, finalmente, a resolução directa do problema 
com forças de massa diferentes de zero, e não só 
a do problema pertubador, correspondente a for- 
ças de massa nulas. 


No problema cuja resolução se apresentou, 
ambos os métodos puderam ser aplicados e a con- 
cordância dos resultados justifica que em ambos 
se deposite confiança. Parece legítimo estender 
essa confiança a casos semelhantes em que só o 
dos elementos finitos seja aplicável. 
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Fig. 25 — Deslocamentos verticais ao longo do estrato AB 
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Fig. 26 — Deformação do conjunto 


Fig. 27 — Deformação. Pormenor 
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DO MUNDO TECNICO 


UM APARELHO DE RADIOGRAFIA 


INDUSTRIAL. 


POTENTE E MUITO FÁCIL DE MANOBRAR 


A utilização do gás isolante SF 6 na fabricação dos 
aparelhos de radiografia industrial permitiu fabricar 
aparelhos «monoblocos» mais leves, mais fáceis de 
manobrar e mais robustos, 


E bem conhecido que o poder de penetração da 
radiação duma ampola de raios-X aumenta com a tensão 
de aceleração dos electrões emitidos pelo cátodo, 
As qualidades notáveis do gás SF 6 utilizado como 
isolante permitiram construir um aparelho aliando 
a muita alta tensão à potência: o Baltographe tipo 
RE 3920. Funciona a 320 kV/14 mÃ e pode radiografar 
chapas de aço muito espessas segundo as normas mais 
severas de controle. 

Está montado num elevador hidráulico móvel de 
comando eléctrico, com a unidade de refrigeração, 
indispensável a um aparelho desta potência, incorpo- 


rada neste elevador. 


O conjunto é de fácil manipulação e pode ser 
suspenso duma ponte rolante. É destinado, principal- 
mente, ao controle das soldaduras de construções 
pesadas, difíceis de deslocar; a sua mobilidade permite 
o posicionamento rápido nos diferentes locais que se 
desejam radiografar. 


Este aparelho é igualmente usado em muitos outros 


casos de controle industrial, como por exemplo o das 


as fundidas. - : 
peças fundid (Informação: S. A, Usines Barreau) 


BLOCAGEM DO ROTOR SOBRE O EIXO EM MOTORES 
ELÉCTRICOS POR UMA RESINA ANAERÓBIA 


Uma nova tecnica acaba de ser adoptada por uma 
sociedade alemã para a labricação automática de peque- 
nos motores eléctricos: esta técnica faz apelo a uma 
resina anaeróbia que permite realizar uma excelente 
acoplagem do rotor e do seu eixo — montados à pres- 
são — e que elimina os riscos de deterioração dos 
elementos a unir. 

Esta técnica elimina a necessidade de utilização de 
peças de oficina muito precisas, quando se faz a acopla- 
gem automática sob pressão. É possível admitir uma 
gama de 1/100 à 4/100 mm entre o eixo e o rotor,o que 
evita o risco de deformar a bobinagem ou de partir 
o eixo no momento de encaixe. 

Esta técnica de montagem é utilizada para uma gama 
de motores cuja potência vão de 1/400a 1/24 CV. Estes 
pequenos motores de corrente alternada destinam-se 
a ventiladores de radiadores eléctricos, blocos de 
congelação e ventiladores de projectores de filmes. 

Antes de passar à operação de montagem procede-se 
ao alisamento das superfícies. Em seguida coloca-se 
uma pequena quantidade de resina no ponto de união 
eixo-rotor. À resina é repartida sobre toda a superfície 
de contactos das duas peças através da rotação que se 
imprime ao eixo, 

A resina fica inerte ao ar, mas confinadas entre 
duas superfícies metálicas, transforma-se numa massa 
dura resistente aos solventes e capaz de assegurar uma 
sólida fixação entre as duas superfícies. 


— Motor shaft 


— Motor rotor 


Retaining 


Compound 


(Informação: Locrrre) 
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ESTUDO DE VIGAS SUBMETIDAS A CARGAS 
PARABÓLICAS 


RESUMO 


No presente artigo estuda-se a viga simplesmente 
apoiada e a viga encastrada sujeita a carga parabólica, 
nomeadamente no que se refere a esforços transversos, 
momentos flectores, rotações e flechas, Apresentam-se tabelas 
e ábacos para o cálculo daquelas características, 


1 — INTRODUÇÃO 


O número cada vez maior de estruturas com 
grandes vãos sujeitos a cargas parabólicas, faz 
com que estas cargas não devam ser assimiladas 
a triangulares, trapezoidais ou rectangulares, dada 
a economia que poderá resultar do estudo mais 
rigoroso do problema. Acresce ainda o facto do 
conhecimento da flecha máxima poder ser, em 
certos casos, decisivo no dimensionamento da 
viga necessária ao estabelecimento do vão. 

Estes motivos levaram à publicação destas 
notas que procuram sintetizar, para as vigas 
simplesmente apoiadas e para as vigas encastra- 
das, o que de maior interesse existe para quem 
tenha de as dimensionar. 

Considere-se um vão de comprimento 1, sujeito 
a uma carga parabólica de valor máximo p e 
sejam x e y as coordenadas correntes da parábola 
definidora da carga, em esquema ter-se-á: 


(xy) 


A lei de variação da carga parabólica será: 


asi q 
pl) =4p (1 É) (1) 
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SYNOPSIS 


The present paper is a study of beams with either simply 
supported or built-in ends, subjected to parabolic loads. 
Shearing forces, bending moments, slopes and deflections 
are dealt with, and numerical tables and charts are pre- 
sented which allow the calculation of those magnitudes, 


O valor da carga parabólica parcial com- 
preendida entre os pontos da abcissa O e x, 
será dado por: 


Po)=[p (x) dx (2) 


ou, atendendo a (1) e integrando entre os limites 
considerados, vem para a carga parabólica parcial: 


Pt) =Tpi[a-2X] (3) 
3 ] | 


o que dá, para a carga parabólica total (x = 1) 
o valor: 


2 
P AA = 
E ELA (4) 


Interessa determinar, também, a posição da 
resultante da carga parabólica parcial, isto é, a 
abcissa do seu centro de gravidade, que será dada 
pela expressão: 


G (x) = (5) 


Atendendo aos valores de (1), (2) e (3), resulta 
depois de integrar e simplificar : 


O DO) ut emma e (6) 


Esta expressão mostra que a abcissa do centro 
de gravidade é independente da carga máxima p, 
como seria de esperar. 


2 — VIGA SIMPLESMENTE APOIADA 
2.1 — Esforços transversos 


As reacções de apoio duma viga simplesmente 
apoiada sujeita à carga parabólica atrás definida, 
são iguais e têm como valor, atendendo a (4): 


Ra=Rg=-——pl (7) 
O esforço transverso, em qualquer secção de 
abcissa x, será dado por: 

T(x) =R, — PG) (8) 


ou, atendendo a (3) e (7): 


1 
T(x) = EL o 1— 6 


2.2 — Momentos flectores 


O momento flector em qualquer secção de 
abcissa x, será dado por 


Mix) =x.Ra — P(x) | *—S6 | (10) 


ou, atendendo a (3), (5) e (7), resulta depois das 
convenientes simplificações : 


aa 2 4 — 


pá X 
a Fe 


X 


tô sea 
pp == ] 


(11) 


O momento flector é, como se sabe, máximo a 
meio vão e tem por valor: 


M o o (12) 


5 
correspondendo a A do momento flector devido 


o 
à carga uniforme p e a E" do momento flector 


devido à carga triangular isósceles de ordenada 
máxima p. 


2.3 — Rotações e flechas 


A determinação da flecha em qualquer ponto 
da viga, poderá ser feita a partir da expressão 


(13) 


Eo)=— [ f des É ai 


em que E e I representam o módulo de elastici- 
dade e o momento de inércia da secção da viga 
de abcissa x, que no caso que se irá tratar se 
supõem constantes em todas as secções. 

O valor do primeiro integral indicado na 
expressão (13) representa o ângulo, em radianos, 
que a tangente à deformada da viga no ponto de 
abcissa x faz com a posição inicial da mesma 
(horizontal na maior parte dos casos). Assim, 
atendendo a (10) e (13) ter-se-á: 


3 4 - 


2 
- 2h 2 * | 
=" [| -2>+> |dx (1 
? (x) pr E x (14) 


Pr x x xé 
og=E | —->+>+kK | (15) 
2 SJ 


em que o valor da constante de integração K, se 
obtém pela condição de anulação de v (x) no 


| 


ponto x == a À donde : 

K, = — e (16) 
e portanto, substituindo em (15): 
pia EL 2 5X 45% a | 09) 
30El / 1 E 


A integração de (17) conduz à expressão de 
flecha, donde, por ser nulo o seu valor na origem: 


Et) = PL 3X EE a e 
“o gp a 


Da expressão (18) se obtém o valor máximo 
de flecha ocorrente na secção média da viga: 


(A ) — pl 
2 5760 El 
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(19) 


2.4 — Tabelas e ábacos 


Os valores das características estudadas poderão 
ser tabelados e traduzidos em ábacos assimilando 
as expressões que as definem ao produto duma 
constante C, dependente só das características 
geométricas e elásticas da viga e dum factor K 
que é só função da posição relativa da secção a 
estudar. Assim, as fórmulas (3) da carga para- 
bólica parcial, (6) da posição do centro de gravi- 
dade desta mesma carga, (9) do esforço transverso, 
(11) do momento flector, (17) da rotação e (18) 
da flecha, poderão ser transformadas conforme 
se indica no quadro I, pela introdução da variável 


X 
k=5. 


No Quadro II, inscreveram-se, para cada carac- 
terística estudada, os valores do factor K e os 
seus acréscimos AK para os valores da variável 
k diferindo de 0,1. 

O ábaco 1 traduz gráficamente a variação do 
factor K relativo a cada uma das características 
que se vêm estudando. 


3 — VIGA ENCASTRADA 
3.1 — Esforços transversos 


Dada a simetria da viga e da carga, os esforços 
transversos são iguais aos da viga simplesmente 
apoiada, donde, devido a (9): 


- à 3 — 
1 x x 
ada | nin + 4 É | (20) 


3,2 — Momentos flectores 


O momento flector genérico é dado pela soma 
algébrica do momento flector isostático (11) com 
o momento de encastramento M,, dado pela 2.º 
expressão (25) abaixo deduzida, pelo que se 
poderá escrever: 


4 


3 
1 x x 


5 | | 


+ E. (21) 
/ 


O máximo momento flector positivo verifica-se 
a meio vão, e tem por valor: 


QUADRO I 
| l 
s | Decomposição 
Características E o | 
O Constante Factor 
K 
Carga parabólica parcial (3) ds | ak —-2k 
MEI Es SEREI DR canas 
Posição do centro de gravidade (6) L | 4 . o : ; 
Esforço transverso (9) 5 pl 1-6K +4k 
Momento flector (11) 9 P a 3(k—2zk + k) 
Rotação E (17) pl. ak — skt , k2 
é 30 EI cs ARE TE Rm À 
Flecha (18) pI ks SE = 
90 El = = 
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QUADRO TI 


Viga simplesmente apoiada 


1 1 1 2 pl p* 
a P 1 2! 3 PI 9 PI 30 EI 90 El 
d Carga Centro de Esforço Momento 
A parabólica gravidade transverso flector Rotação Flecha 
E) P(k) G(k) T(k) M(k) P(k) E(k) 
O | ju meme E RC SR (E ia E 
> K AK | K (AK K | AK K | AK REL k !1 AK 
() O ps Õ — | +1,000 |—— 0 — 1,000 — O ai 
fu À 2088 ==) qua AS | — 0,056 +0,294 | +0,050 |——— | 40,295 
0,1 |+0,028|——— | +0,182 + 0,944 | — -+0,294 |——| — 0,950 |— | +0,295 |— 
——— | — — R — +=0,130 | — — 0,152 |—— | 40,264 |- +-(,143 0,266 
0,2 | +-0,104 |—— | +-0,262 +0,792 |——— | 4-0,558 — (1,807 — +0,561 — 
— | — +-0,112 |———| +0,126 |-————| — 0,224 |——— | +-.0,204 ——| +0,221 —| +0,211 
0,3 | +0,216 + 0,388 + 0,568 |-— | 40,762 |— — (1,586 |— 0079 | 
—| 0,186 |————| +0,121 —| (0,272 |— +0,132 |— +0,218 | — +0,135 
0,4 | +0,852 —| -+0,509 |-— +0,296 |—— | 40,894 —| —0,308 cmi DOT | 
pe -++-0,148 +(0,116 |— —[(,296 |——— | +.0,0492 +0,308 |—— | 41.0,046 
0,5 | +0,500 |— 40,625 | — O a + (1,936 |—— — +0,953 
Es —| +0,148 |————| +-0,108 = —(),296 |— — (3,042 -+0,308 |— | — 0,046 
0,6 | +0,648 | | +0,783 — 0,296 —| —+(,894 +0,808 —| +0,907 | 
Ei | +0,136 —| 0,098 —(1,972 |— —| —0,182 | +0,278 | -—| — (1,135 
0,7 | +0,784 | + 0,881 — 0,568 |— — (0,762 (1,586 |— o e É ÉS pe 
—— | «++(,112 |— + 0,083 |—— | —0,224 — (1,204 |— — | +0,221 |——| — 0,21 
0,8 | +0,896 +0,914 — 0,799 —| +0,558 10,807 || Aos E 
+0,076 |— —| +-0,061 |— —0,152|———— | —0,264 |—— | +0148|-— | —0,266 
0,9 |+0972 | +44097,5 |— —(,944 —— | —Q,29M4 || 40,950 |———— | +0,995 | 
— | +0,028 |-— ++0,025 —| —(,056 —| —0,294 |— | +0,050 |—————| — 0,295 
= 1,000 | es me O (O) a 


1,0 | +1,000 |— —+1,000 | 


l d 
M(—)=——pEÉÊ. (22) 
2 80 


Como verificação constata-se que a diferença 
dos momentos flectores a meio vão, dados por 
(12) e (22) é igual ao momento de encastra- 
mento dado por (25). 


3.3 — Rotações e flechas 


De acordo com a expressão (13) e com as 
restrições indicadas em 2.3, a rotação genérica é 


dada por: 
1 É E x 
o (x) =— E oi 
| (x) El / | 3 l + 
x* 4 
+ 1º +M,. | dx (23) 
284 


cm) 41,000) [mm 


em que M, é o momento de encastramento 
integrando (23) vem: 


E 3 2 4 
Go) mad (= BE. 
EL | 3 à2f AM 
És - (24) 
" A fd 


devendo verificar-se 7 (0) = 0 e v ( ) = O 
2 
por condições de compatibilidade. Daqui resulta : 


2 
K, ='0 M=— E (25) 


Substituindo (25) em (24) vem finalmente. 
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ÁBACO 1 — Viga simplesmente apoiada 


RE E 
Cons siansasnasdiaaa spdains as IPRHAS RR 
10 qro ONES NANESES FETO Leer 
E O a 
ERA NEH f ps A a 
cu EorsnsssEadE SS ESEsI SSESS NNE A 
9 ENE ER AR dipastiaá Fam sida, CORsMEaMaRaSaSs qÁES feER asas tosa: 
CORRE SOS S Sopa Ea qa io PRGREEREEECE EEE ONG RE 
ssssnseEs ABRA PI EE NNE 
sous es ESET idgrai Pe htsdisdi GENDER S TANTA ame. Pr 
oops Ni A ni oa srs 
Cqssasssos RN Hs EA RNP EIS 
COSES Ce Sos SESSEESSSSER QUES /AESTSSESES SRS REER ESSES pesas é cpas' fasssesispsssmnsansõass 
EEERParEtREccodbadaa Ms ASS patos os ssdEcssSecsiEssso oaiso ÁSIE, Leds cretpastnicmeeteateas: 

[6] 7 CORRE SSGOS EnSaA SSB SERNADAD sta na td teres sbemaaman aaa anos 


EE / rendas sonsas CDS NnEa" daunsepaas VaUDa Posso Udo Sme ue Sana 
ESPEP CINE 
ie OR CT ARE O De o tm o q 1 pd e 


j 
| 
: 
E 
no 


o eus AE 

did EEssEsEsssssssssassoo Srrtssa cases Eorastoasosraforeoasr Snstaspadtaas jstsstea, StEEsiscasceaansesss 

COST SrESARARaEaSS AEssEatE SEnEMEaERE EA io S87 (POSSR SARA VERAS PosES ceces secas 

ERRREERERtEatESar HOotoCroai rEStEdiss dtaaiEa DS nsaicaroa too Aotina pesto, ospostcacsscáseas 

SEE SERA SIDE 

EEIAESEEiSaso ObpoEsessptase dos pstadisS (é PEaGisadsd atodt cá poa Toa RE SESS ses tismsamiaiadass 

Faso tECEC EEE Codpesiscestadios, cafssfatas! jo fpadi sa ná pads MfEss SEsstastRatEat, HESSE SRSSEStES 

O 4 sesas innso DA GUaa Rana RESRN E MERO VON DNA dENSSDANANA RENNA DOR DUNNUS RUNAS RannRnDa ara Sana RSannanNs 

. DO OD O 4 A eme aaa 42.e s ra ad HE | RO À SA Sa RES ns 

ERR DES ita fecoEONEAEEE EEcE CERCAS ERC CHER ERR 

FREE EEE iss Resid Ea nar Sbt root tEaissssdn testestsatcss 

o,a EEE E FS Ed Pass saRRA LAFER fis saaaasRREaEEasifsssssa ara 

Es A ELE dEnd setas tan Eapádadss cotsdnapadoassrosERSsSESRStEGAS canstasss 

Pavanas nana” EEN TUIS CESSSNNES IE A ESSES MASEI ESSES SasES ssa =s casa (à coca acess 

RO RAE ORAR ERREEEE sa cce 

0,2 EEE RR Rs 

EEE psi AA ERA 

a CESSA SNS DERNNRRNHA FRGGEEACEE II eai bh 

a 

ERP AA EE 

É decano RaaR dora tasda dios toctos test tas toscadlastcstas tos 

o Ee E E ML a 
CEsSaE EH 0,2 E HH 0,4 H EH Hr 0,8 HH 0,9 Es 1,0 (k) 

EEE H SE TITI : DT mtH , mt = Hi entttttrerottrto 

Gis ii na E 

OM faso po NE pe seaaa entao pn 

hd [] n numas CEU anansnREnanannaNanaa 

Esp a o SS 

EDEMA - sa: od E + DO E na a a a 

| pa da Ri a 

EFA Doiassasass anda asEasnanas 

E aee qnto una Ae 

| HE it; ua 8 EE Rondonia 

03H Ri] E : SRSPIERSES TENS: Peas pacas de sas cenas 

en o 

fersses E ts E ui AopraniEsiiasstasEss 

a BR ns 

AR flo Haia RE 

GUERASE 7 HH 

EESss Bi MN ptsstpasistiascadsdteass 

| EEE Esses Ei EEE do) fado ESaasEsigacRads CERSs Essas ces 

RA Ro 

anne Aos 

ddsas sosamss f à HH no: Esnasuunas 

EoSESNSSES Gases nasas! dessas se à Qfasisddos issssisdos nana canas 

nan at A 

disnenssss isdan nas dbesasia HEHE | Rosi spasa setas ineo 

EoRSEsiaEs Pers PESR ds Eatasas E EEREENIO Ao OR GER 

fas ARA ERaRaAANNRRA VRRSaSGS SS cosas 

Eostecarcsanosiscatioss GttnatisaEs ditEsa tidas ioat nah Sbt catrrasicesttesi£os 

ERES RR E 

EoiEoPosao tes dotiatass d pisos cilaiaisaiaicascsisicatsstacasiadadictas RN HERE EEE 

oops fa] — convENÇÕES Nou 

HH HT a ISERSES AS asanssas 

RSI io | canoa PARABÓLICA A 

ERA HH] CENTRO DE GRAVIDADE — — — | ESION E erre 

CSSSSESSS Ao 

-0,9 Eat] ESFORÇO TRANSVERSO — —— [pe 

paca ngicifi o] MOMENTO FLECTOR —.— | NE 

CLA ERsddot] ROTAÇÃO — e ER EEE 

A EL E ERES FLECHA EHEHE 

1,0 SUNSRARNNARCRNanaaNaaNaaNGNaSS DONS NSSAERRASARnanaanaaaNaano 

EORtEA TERES rESPoSioatoSto PES fassss pares posPestosapaçaspessosnoasaspasnas ta diadfiatasistastiiatoss 


Pi : x O valor máximo da flecha obtém-se a meio da 
g=E. |2Ê-sÃos 
30 EI Ú 1 viga e tem por valor: 
oO (26) 
posa 
Ê a | 13 pI 
(sm 
- 2 5760 El 


A integração de (26) conduz à expressão da 
flecha, donde, por ser nulo o seu valor na origem: 


- 3 
file pl. à , =. 3.4 — Tabelas e ábacos 
90EI 2 / 
5 65 (27) 
Alicia, De maneira idêntica à indicada em 2.4 orga- 
P ds nizaram-se os Quadros Ill e IV e o ábaco 2. 
QUADRO III 
| cá Decomposição 
| 3 | —— 
Caracteristicas E 
| o Constante | Factor 
Co | K 
Carga parabólica parcial (3) - pl sk — 2lé 
Posição do centro de gravidade (6) : l 4k — 3k 
- | 3 — 2k 
a | E | 
Esforço transverso (20) ; pl 1-6k + 4k 
Momento flector (21) , pl — T+ Sk — 10k' + 5k* 
e se | 
- | P o 5 4 o 
Rotação (26) ES =) Ss(2k — 5k + 5k — 2k) 
OD e ERES E Rca sas o 
p | 2 3 5 6 
Flecha | (27) 450 El | S(3k — 5Sk + 3k —k) 
| 
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QUADRO IV 


Viga encastrada 


Cons- 2 1 1 5 |” Jú 
Zz É ic A 9 E. Pp 
Ea 3 Pl 2 | 3 P1 3 PI 150 EI 450 EI 
, Carga Centro de Esforço Momento R lech 
o parabólica gravidade transverso flector ração FRA 
; | Coma Gt) T(k) Mík) EU) Eq) 
Õ E = a a 
> K [AK | K AK K | 4K K SK K A K K A K 
0 O [loss q — 41,000 |—— | —1,000 =D “A E, 
eme + 0,028 TO 152 ———, —0,056 |——| +0,490 |— | —0,752 |—— 1 40,125 
0,1 | +0,028 — — | +0,182 + 0,944 —U,510)/— gs |— 140125! 
Do | +0,076 |—| +0180|— | —o 152 |— | +0,438 — 0,285 |— | 40,979 
0,2 |+0 104 |— | +0,262 |——— | 40792) | -0,079 — 1,087 + 0,404 
Do rr) FOZ |———| +0,106 |———) —0,294 |—— — | rosas |— +0,109 |-— >! 0,804 
0,8 | +0,216 |——— | +u,888 |— +0,668 |— | +0,270]|—— | 0928 |—— | +o708|-— 
—>—|— ——| +0,136 |—— — | +0,1921 |— —0,212 |—— +0,218 |—!| 40,890 |—— | 40,295 
0,4 |+0,3552 -———| +0,609|— | +0,296 | | ioass|— | osss|- +0,988 |-— 
—— | | +0,148 +0,116 |———| —0,296 |———.| +0,074 |— | +0,588 |-———— +0,088 
0,5 | +0,500 |—— | +1.0,625 |— O |-——l+o56 | q +1,016 |-— 
pn] TOt48 |——| +0,108 |——! 0,296 |— — | —o07 | ++0,588 |— | — 0,083 
0,6 | +0,648 |———)| +0,783 |-—— 1 0,996 |—— | 40,488 |— -+0,538 +0,988 
o [07 | +0,186 | +0,098 |-— | 0,972 |——— | —0,218 +0,890 — 0,225 
0,7 | +0,784 |— | 10881 —0,568 |—— | +o,270 |— — | 40,998 0,708 
ad RR BET TT —0,224 |— | —0,8492 +0,1 —[(),804 
-+0,896 |—— 1 40,914 — 0,192 — 0,072 -+1,037 0408 | 
+0,076 1|— | — 0,152 —0,438 |—— | 0,285 |— — 0,279 
+98 | +05, — 0,944 — 0,510 +0,752 +0,145 
rn] 0,028 |— — +0,025 |—— | —0,056 |— | — 0,490 | — —(),752 —| — 0,125 
1,0 | 1,000 | 1,000 | — 1,000 “100 o | Do 
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ÁBACO 2 — Viga encastrada 
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] e e 
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C. D. U. 536.25:532.542.4 


DeLaano Domingos 


Efeito da variação axial da temperatura ou do fluxo 
térmico no coeficiente de transmissão do calor por 
convecção forçada num tubo 


Técnica No. 385 — XLIV — 2. 4969, pág. 225-231 


Discute-se a importância das condições na superficie de 
permuta — temperatura ou fluxo térmico imposto — no 
valor do coeficiente de transmissão do calor por convec- 
ção, num tubo, os quais podem diferir substancialmente. 
Utilizando propriedades das equações fundamentais mos- 
tra-se como se podem calcular os coeficientes quando as 
condições impostas variam axialmente desde que seja 
conhecido o seu valor (por via teórica ou experimental) 
para condições uniformes. 


Aral ra Ma aa a ia a ea e hs 


CarLos PonrTELA 


C. D. U. 681.14 :5148,3 


Calculadores analógicos — aplicações a sistemas 
mecânicos 


Técnica No. 385 — XLIV — 2. 1969, pág. 233-261. 


Este artigo trata da aplicação de calculadores analógicos 
a sistemas mecânicos. São analisados os seguintes 
aspectos : 

— Análogos eléctricos. 

-— Tipos de análogos eléctricos. 

— Circuitos lineares de parâmetros concentrados, 

— Análogo da relação entre forças e deslocamentos 
generalizados em elementos actuados em dois pontos. 
Aplicação a sistemas complexos. 

— Aspectos conceptuais e metodológicos do cálculo 
analógico. 

— Análogos de alguns sistemas mecânicos simples em 
regime estático e dinâmico, 

— Análogos de uma viga em regime estático e dinâmico. 

— Análogo de uma barra sujeita à torção em regime 
dinâmico. 

— Regimes livres. Frequências próprias. Representação 
modal, Análise de vibrações. 

— Análogos de elementos não lineares. 

— Análogos de placas em regime estático. 

— Análogos aplicáveis a problemas de condução de calor. 
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Resumo dos artigos publixados na Técnica n.º 385 


Ano XLIV — Fevereiro 1969 


COCEIRA E CEC CU... 


C. D. U. 539.31:519.34 
519,34 
551.434 


E. R, ve Arantes E OLIVEIRA 


Análise elástica de vales e taludes | 
Técnica No. 385 — XLIV — 2. 1969, pág. 263-279. 


Assimila-se a crusta terreste a um semi-espaço elástico 
onde se instalam estados planos de deformação. Formula- 
-Se assim um problema de elasticidade plana relativo a 
um semi-plano indefinido, 

Considera-se a resolução deste problema por dois méto- 
dos numéricos, o método dos elementos finitos e um 
método Integral, e comparam-se os resultados obtidos para 
um caso concreto, o da determinação das tensões num 
estrato argiloso que ocorre no talude das Portas do Sol 
em Santarém. 
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C. D. U. 624.042.2 


Penes Roprraves 
Estudo de vigas submetidas a cargas parabólicas 


Técnica No. 585 — XLIV — 2. 1969, pág. 281-288. 


No presente artigo estuda-se a viga simplesmente apoia- 
da e a viga encastrada sujeita a carga parabólica, nomea- 
damente no que se refere a esforços transversos, 
momentos flectores, rotações e flechas. Apresentam-se 
tabelas e áhbacos para o cálculo daquelas características. 


.« 2.0. 


7,6 AA OA 


A 
CO ra a e ana aaa... 


289 


| 
aaa 


Synopsis of articles publ 


ce s0.00. 


; XLIV — Febr 
OO AAISEKEKKSESA NESSES AAKAÃ secrtriTerE Te as 
$ 


“eee sanu.. “once aca non nn asas as 


WMM: 
 —-= .. 


ao 


“qu d = Dt pe cm dp std tm al) a pal at A ima ab es mm So Se a el PS Rd 


= 


10400400 + 


E. R. pe Arantes E OLIVEIRA 
519.34 
551.434 


Elastic analysis of valleys and slopes 


Técnica No. 385 — XLIV — 2. 14969, pp. 263-279. 


By identifying the earth crust with an elastic half.space 
subjected to a plane state of strain, the elastic analysis 
of valleys and slopes can be converted into a plane 


elasticity problem defined on a half-plane. 


The author discusses the solution of this problem by 
two numerical methods, the finite element method and 


an integral method, and compares the results vielded 


by both when applied to an actual case : the analysis of 
stresses in a clay layer which occurs at the Portas do 


Sol slop in Santarém. 


U. D. G. 624.042.2 


Peres RopriGuEs 
Study of beams subjected to parabolic loads 
Técnica No. 385 — XLIV — 2. 1969, pp. 2814-288. 


The present paper is a study of beams with either simply 
supported or built-in ends, subjected to parabolic loads. 
Shearing forces, bending moments, slopes and deflections 
are deglt with, and numerical tables and charts are pre- 
sented which allow the calculation of those magnitudes. 


290 


U. D. C. 539.31: 519,34 


PUT PPT TT EV e 2"0-0=0 erre + ——— e... 


e e O Tee mem > e 
COCCOCCCLCCLCLCCLCLLCCLLLLCCLLCOCLCLCEr E Cena ra a a e na E Rn a sn 4 


e 


“.«.. 


TFPTTTTTTTTTITTATTALAAAAKKKKKKKKAAKKAMKKMKAKEKS KA EE O OO O O CR 


.......... e... 


hed in «Técnica» nº 385 


uary 1969 


O. D. C. 536.25 :532,54244 


Dec.gano Domingos 


Effect of the axial temperature variation or of the 
thermal flux on the heat transfer coeficient by 
forced convection in a tube 


Técnica No. 385 — XLIV — 2, 1969, pp. 225-234. 


The importance of boundary conditions — uniform tempe- 
rature or heat tlux — on the heat transfer coelicient by 
forced convection in a tube is discussed. 

Starting from the properties ofthe fundamental equations 
for heat transfer it is shown how the heat transfer 
coeficient can be found for axially varsing boundary 
conditions when its value js known (experimentally or 
analyticaliyv) for uniform temperature or heat flux. 
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U. D. C. 681.14: 518.3 


Analogic computers — Applications to mecha- 
nical systems 


Técnica No. 385 — XLIV — 2, 496º, pp. 233-264. 


This paper deals with analogic computers and their 
application to mechanical systems. The following aspects 
are presented ; 

— Electrical analogys. 

— Types of electrical analogs, 

— Linear networks of lumped parameters. 

— Analog of the relation between generalized forces 
and disvlacements in elements actuated in two points. 
r pplication to complex systems, 

— Conceptual and methodolosic aspects of analogic 
computation. 

-- Analogs of some simple mechanical systems in static 
and dynamic regime, 

— Analogs of a beam in static and dynamic regime. 

— Anelogs of a beam submited to torsion in dynamic 
regime. 

— Free regimes, Natural frequencies. Modal represen- 
tation. Vibration analysis. 

— Analogs ot non linear elements. 

Analogs of plates in static regime. 

— Analogs applicable to heat conduction problems. 
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ARTIGOS DE REVISTAS 


C. D. U. 537.77 :541,1 — 145.1 


The standard free energies of solutions of anhydrous 
salta in water — D. 4. Johnson. 

Journal of Chemical Education, 4-68, vol. 45, págs. 
236-40, t. 4. 


Estudo sobre o tratamento das variações de energia 
livre com o tamanho do catião, com predição de várias 
tendências relacionadas com este tratamento. Faz-se 
também um breve estudo relativamente ao tamanho 
do anião. 


C. D. U. 539.16 — 544.62 — 545.844 


Non-radioactive physical methods of analysis dvring 
gaseous reprocessing of nuclear fuel, — W Goossens. 

Industrie Chimique Belge, 7 e 8-968, vol. 7-8, págs. 
643-654, t. 33- 


O control «on-line» de misturas gasosas durante a 
extracção gasosa de combustível nuclear é estudado 
para um «flow-sheet» de CIF;. 

Os métodos físicos não radioactivos que podem ser 
utilizados para analisar as misturas de gases corrosivos 
em causa são estudados segundo o seu aspecto de con- 
trol contínuo ou descontínuo, 

Os resultados recentes permitem constatar que a 
espectometria de infravermelhos e a cromatografia 
gás-líquido são os mais interessantes, 


C. D. U. 541.124/.128 — 145 


Some aspects of chemical dynamics in solution — Jack 
Halpern. 

Journal of Chemical Education, 7-968, vol. 45, págs. 
373-80, t. 6. 


Estudo das diferenças fundamentais de dinâmica 
química em relação à fase gasosa, Estudo da solvatação 
dos reagentes e / ou das espécies intermediárias in- 
cluindo o complexo activado, e participação do solvente 
na reacção, como agente nucleófilo. 


CG. D. 0. 511.127.1.04 
Chemical reaction cross sections and rate constants 
— E, F. Greene and A. Kupermann. 
Journal of Chemical Education, 6-968, vol. 45, págs. 
361-g, t. 6. 


Estudos de velocida desde reacção e estrutura mole- 
cular em reacções cruzadas bimoleculares, em fase 


gasosa. Utilizam-se experiências de colisão bimolecular 
num gás, enunciam-se métodos experimentais para a 
determinação da secção recta de reacções e aplicações 
destas reacções aos líquidos. 


C. D. U. 541.128:532,7 
Catálise y cinética química — Anthony Tucker. 
Ion, 10-968, vol. 28, págs. 712-4, t. 327. 


Estudo das enzimas e das suas propriedades como 
catalizadores e aplicações industriais das mesmas* 
Acção dos catalizadores durante uma reacção, e estudo 
do mecanismo pelo qual estas substâncias aceleram a 
reacção. 

C. D. U. 545.242 

Mechanism of oxidation — reduction reactions — Flenry 
Tavbe. 

Journal of Chemical Education, 6-968, vol. 45, t. 7, 


pág. 45361. 


Considerações sobre número de oxidação e estados 
de coordenação, estrutura electrónica e estado de 
coordenação, complexos activados por transferência 
electrónica, estrutura eléctrónica e mecanismos e 
transferência electrónica por ataque remoto. 


G. D. U. 514.64/65:515 


Correccion de influências mútuas elementales en mues- 
tras naturales eplicando corficientes matemáticos de 
correccion para determinaciones quantitativas por 
fluorescencia de rayos x — /. L, Sagrera. 

lon, 12-968, vol. 28, t. 329, pág. 856-61. 


Aplicação das diversas fórmulas matemáticas de 
convecção de influências mútuas elementares em 
técnicas quantitativas por fluorescência de raios x, 
dando uns exemplos de aplicação que permitem 
comparar a validade dos ditos coeficientes matemáticos 
de convecção. ; 
CG. D. U. 547.971 
The methanolysis of substituted methyl benzostas — 
J. T. Glover and C. W. Peterson. 

Journal of Chemical Education, 4-968, vol. 45. t. 4, 
pág. 241-43. 


Trata-se do estudo de uma reacção de substituição 
ilustrando as relações de variação linear da energia 
livre do processo. Esta experiência é feita com metil 
benzoatos substituídos em presença de álcool metílico 
em que o primeiro tem um carbono radioactivo que 
nos permite seguir a reacção. 
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TÉCNICA — XXIV 


C. D. U. 551,482.212.3 


Avaliação do caudal sólido em escoamentos unidirec- 
cionais — Z. Veiga da Cunha. 
Fomento, 1968, vol. 6, pág. 21. 


Nesta publicação passa-se em revista a aparelhagem 
que tem sido utilizada na medição directa ou indirecta 
de caudal sólido por arrastamento e em suspensão — 
fazendo-se a apreciação de cada tipo de aparelho, 
Descrevem-se técnicas de medição de caudal sólido, 
analisam-se problemas do número e localização dos 
pontos de recolha de amostras e ainda da frequência 
de recolha em cada ponto, 


C. D. U. 552.1 


Critêre de rupture des roches fragilles — V. S. Parate. 


Annales. 1-969, vol. 253, pág. 149. 


O comportamento da rocha é dúctil, quando subme- 
tidas a elevadas compressões médias e frágil nas con- 
dições vizinhas das atmosféricas. Na maior parte das 
aplicações, o engenheiro tem que enfrentar problemas 
de rotura frágil ou porque deseja evitar a rotura ou 
porque deseja provocá-la. Estabelece-se a distinção 
entre dois critérios: 

— quando, pelo menos, um esforço principal é uma 
tracção, a rocha fractura por extensão e obedece aum 
critério de tracção máxima; 

— quando os três esforços principais são complexos, 
a rocha obedece a um critério de corte máximo. 


C. D. U. 621.165-52 


Mise en service d'un dispositif de commande aunto- 
matique de groupe turbo-alternateur au Canada — L ger 
Dieter. 


Rev. Brown Boveri, t. 54 (1967), n.º1, p. 3-8. 


O autor descreve as características e a constituíção 
do dispositivo electrónico turbomatic e indica as fases 
da sua montagem. 


C. D. U. 621.5.052.63 : 621.398 


Appareils à conrants porteurs à deux bandes latérales 
ZSB 400 pour télóphónie et téléopération sur les lignes 
à haute tension — Aigner, W.; Domes, R. 

Rev. Siemens 27 (1969) n.º 1, pages 24-30. 


O artigo descreve a constituição e o funcionamento 
do aparelho e indica as suas aplicações, 


C. D. U. 621.311,22-52 


Le Decontic, un systême de commande décentralisé pour 
Vautomatisation des centrales — Môtz Friedrich 
Rev. Brown Boveri, t. 55 (1968), n.º 1/2, p. 25-31. 


Após uma descrição sucinta do sistema de comando 
o autor descreve o funcionamento dos módulos 
utilizados neste sistema, assim como a sua montagem, 


C. D. U. 621.313:333.620.179 


Estatística des correntes em vazio de motores assin- 
cronos— Franklin Guerra, 
Electricidade 1-2-69g n.º 57 pág. 47-48. 


O autor faz uma análise estatística das flutuações 
em relação à média, que se verificam nos ensaios de 
verificação, à saída duma linha de montagem. 


G. D. 0. 621.313:562 
Une nouvelle serie de moteurs shunts triphasés à 
collecteur — Blunsch: Adolf. 
Rev. Brown Boveri, t. 55 (1968), n.º 1/2, p. 20-25. 


O autor descreve a nova gama de motores. 


C. D. U. 621.317.79:621.3,518,46 
Dispositif de mesure du glissement, utilisant des cir- 
cuits intégrés — Miiller Edmund et Skuin Petar, 
Rev. Brown Boveri, t. 54 (1967), n.º 4, p. 187-189. 


Neste artigo descreve-se um aparelho que, por meio 
de circuitos integrados, a partir da frequência do 
«stator» e da velocidade de rotação do rotor, elabora 
um trem de impulsos representando a diferença de 
frequências e por conseguinte, a frequência do escor- 
regamento. 


C. D. U. 621.316.364 


Multi Compact—un noavean systême intégré d'appareils 
de connexion et d'équipements de connexion — Paath, S. 
ASEA-Revue 40 (1968) :6, p. 123-125. 


Graças a uma coordenação lógica, a sua adaptação 
sistemática ao meio e a uma execução compacta, o sis- 
tema «Multi-Compact» permite fazer uma instalação 
com menos custo. 


GC. D. U. 621.316.1-523.8:621,313.362 
Régulateur à trois paliers pour moteurs shunt triphasés 
— Neumann, B.; Waldinger, HH. 
Rev. Siemens 27 (1969) n.º 1, pag. 30-33. 


O artigo descreve o aparelho, indica o modo de 
funcionamento e as suas passibilidades de emprego. 


G. D. U. 621.38-112 
L'insensibilité aux perturbations des systemes électro- 
niques à éléments de construction — Morgenthaler Max 
Jakob, 
Rev. Brown Boveri, t. 54 (1967), n.º 4, p. 1156-161. 


No presente artigo analisam-se os factores essen- 
ciais da insensibilidade às perturbações de sistemas 
lógicos electrónicos à base de elementos de construção, 
referindo-se particularmente ao caso de circuitos inte- 
grados. 


C. D. U 621.99:621.317:681,3 

Equipements automatiques de contrôle pour la techni- 
que des télécommaunications — Kemp, W. R. 
Rev. Siemens 27 (10969) n.º 1, pág. 17-20. 

O artigo trata da construção de aparelhos automá- 


ticos de controle e de medida. Descreve entre outros 
o conjunto automático de medida «PEGAMAT», 
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TECNICA — XXV 


C. D. U. 621.395.2 


Grandes installaticns tvléphoniques à postes supplé- 
mentaires ESK S(00 E réalisées en technique Siemens 
Crosspoint — Aurz, W., 

Rev. Siemens 27 (1969) n.º 1, pag. 21-24. 


O artigo descreve as características fundamentais 
das instalações ESK 3000 E: Construção e constituição 
do sistema, rede de conexão e comando electrónico. 


C. D. DU. 621.396.7:620.19:550.3 


Extensicn de la radio-station terrestre de Raisting pour 
le trafic des tólécommunicesticns pzr la satellite Intel- 
sat UI — Graf, H ; Schiirkimper, K -H. 

Rev. Siemens 26 (1968) n.º 5, pág. 147-152. 


O artigo descreve a trunsformação da estação de 
rádio Raisting | para exploração por intermédio do 
satélite Intelsat Ill e faz referência à montagem da 
antena da estação Raisting IL. 


C. D. U. 621.742 + 621.744: 620.1 


Novos aparelhos de control de areias expostos na GIFA 
1968 — O. Hikl. 
Revista A. P. F., n.º 31, 968, pág. 21. 


O artigo descreve uma nova gama de aparelhos de 
control de areias, cuja concepção tem em vista obter 
resultados laboratoriais adaptados à realidade das ope- 
rações na moldação solucionando algumas contradições 
que se verificam coin os ensaios feitos nos tipos conven- 
cionais de instrumentos. 


C. D. U. 621.743 + 621.744 


Máquinas-ferramentas modernas na indústria de fundi- 
ção — /. Fallows. 
Revista A. P. F., n.º 31, 968, pág. 9. 


O autor dividiu esta comunicação ao 2.º Congresso 
Nacional de Fundição em duas partes — na 1.º relata a 
razão do atraso da mecanização no sector da indústria 
fundidora e apresenta um caso concreto de uma colabo- 
ração bem sucedida entre um fabricante de máquinas 
para a fundição e uma empresa fundidora; na 2.º parte 
descreve a máquina de moldar sem caixas, denominada 
Disamatic. 


C. D. U. 62407:001.57 +518.5 


Utilisation du calcul numérique et des modéles réduits 
pour l'étude de constructions — £. Absi e Ch. Bonvalet. 
Annales 2-69 vol. 254 pág. 259. 


O estudo das estruturas pode, em casos difíceis 
recorrer a duas técnicas complementares: 

1.º — O cálculo electrónico 

2.º — ()s modelos reduzidos 

O primeiro autor desenvolve diferentes aspectos e 
aplicações do cálculo electrónico. O segundo autor 
apresenta um certo número de exemplos de estudo em 
modelos reduzidos. 


C. D. U. 24.07:519,5> 


Metodo: matriciales para mecanizar el calculo de estru- 
turas — Frederico Alicart. 
Revista de Obras Publicas 1-69 vol. 3545 pág. 3. 


Expõe-se em síntese, o processo de mecanizar o 
estudo de uma estrutura, com base no método da rigi- 
dez para o que se requerem conhecimentos sobre 
programação; se as estruturas excedem a capacidade 
de memória de um computador, podem ser estudados, 
decompondo-os em subestruturas elementares, 


CG. D.U. 624.012,45 


Hormigon armado: abacos para el calculo de secciones 
rectangulares sometidas a compression compuesta y a 
flexion compuesta — /racísco Ruiz Garcia. 

Revista de Obras Públicas 23-69 vol. 3046-3047 
pág. 137. 


O autor encontrou umas relações lineares entre os 
espaços a que estão submetidas as secções e as 
características geométricas e de tensão das mesmas. 
Estas relações permitem a construção de ábacos de 
fácil manejo e extrapolação, que permitem determi- 
nações imediatas, 


CG. D. U. 624,012.45:624 026: 513.4 


Escadas helicoidais e para-helicoidais de betão armado; 
método de cálculo gráfico-analítico— /. A. Ribeiro 
Sarmento, 

Engenharia (F. E. U. P.), 1-2-968, pág. 5. 


O método que se descreve é aplicável a escadas 
de apoios perfeitamente encastrados, apresentando, 
em projecção horizontal, simetria relativamente ao 
eixo médio da linha dos apoios. 


C. D. DU. 624.13851:620,1 


Alguns resultados de ensaios de laboratório de solo- 
“Cimento com «Muceque» — FI. Novais Ferreira, 
Fomento, 1968, vol. 6, pág. 3. 


Dentro de um programa geral de estudos de solo- 
-cimento com solos de Angola, foi realizada uma. 
primeira série de ensaios com «muceque» — Solo 
argiloso avermelhado — constituída pela obtenção de 
curvas de variação de baridade, resistência à com- 
pressão e CBR, com o teor de humidade para várias 
dosagens e várias energias de compactação, Avalia- 
ram-se as coordenadas dos máximos dessas curvas, 
coordenadas essas que se estudam neste trabalho. 


C. D. O. 663.63 

El problema economico en la depuracion de las aguas 
residuales — Eduardo Dewisme. 

Revista de Obras Públicas, 1-969, vol. 3045, pág. 29. 


O autor pretende dar a conhecer algumas indicações 
sobre o custo de construção e conservação das plantas 
de depuração de águas residuais urbanas, assim como 
os critérios que vêm sendo utilizados noutros países.. 
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CONSTRUTORES 


CALCINA 


E hoje em Portugal o melhor ligante hidráulico 


para argamassas, rebocos e alvenarias. 


Preliam, pois CALCINHA 
a 


Entregas permanentes de dia e de noite 
nas Fábricas de Cimento “Tejo” - Alhan- 
dra e do “Cabo Mondego” - Figueira da 
Foz, e através da vasta rede de reven- 
dedores espalhados por todo o* País 


TÉCNICA XXVII 


SDL. 


ESTUDOS 


Porto 
R. de Sá da Bandeira, 651 - 4 
Tel. 27013 


E INSTALAÇÕES ELÉCTRICAS . 
Lisboa 


L. Campo Pequeno, 21 - 4º Esq. 
Tel. 7680964 


Projectos e instalações 
eléctricas de: 


Centrais 


Subestações 


Postos de Transformação 
Fábricas 
Edifícios públicos 


Blocos residenciais 


1 


Es | A , k 
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CELAS ABERTAS DE ALTA TENSÃO 


Miller 


Níveis e Teodolitos 


Os melhores preços do mercado 
Assistência Técnica Garantida 


Representante exclusivo: 


Nucleon - Equipamentos de Precisão, Lda. 
Av. António Augusto de Aguiar, 165, r/c - Dto. 
Telefone 530105 
LISBOA 


TECNICA XXVIII 


TOPOGRAFIA GERAL 


2 º Volume 


Pelo Eng.º A. C. XEREZ 


Nova edição revista e actualizada 


Preço 150$00 


Desconto : 10º/ aos assinantes 


Pedidos à Técnica 


LISTA CLASSIFICADA DE ANÚNCIOS 


Neste número continua a «Técnica» a sua secção de publicidade 
— a lista classificada de anúncios — destinada a facilitar a procura 


por parte dos clientes. 


Estas páginas amarelas ficam abertas aos nossos leitores, que nelas 
podem anunciar a sua firma e endereço, ao preço de 50$00 por número. 
Quanto aos nossos anunciantes, é com muito gosto que lhes concedemos 


grátis a inscrição na lista. 


CONSTRUÇÃO CIVIL 
CONSTRUTORES CIVIS 


— Construções Técnicas, Lda. 
P. do Município, 13-3.º — Tel. 36 6506 


— Fundações Franki 
R. Braamcamp, 13-3.º — Lisboa — Tel. 34112. 


— Soc. Const. Amadeu Gaudencio, Ld. 
R. Dr. Alexandre Braga, 4-A — Lisboa 
Telef. 43191 — 92 — 93. 

— Soc. de Pré-Fabricação e Obras Gerais 
Novobra 


Av. Estados Unidos da América, 100-5.º D 
— Telef. 77 48 32 / 7664 46 — Lisboa 5 


— SOMAGUE 
R. Carlos Testa, x — Tel. 5624 81 — Lisboa 1 


ESCAVADORAS E TRACTORES 


— Guedes & Almeida 
Rua Aurea, 181-2.º— Tel. 3278 45 — Lisboa 2 


— S, T. E. T. — Soc. Técnica de Equipamentos 
e Tractores, Ld.' 
Prior Velho — Sacavém — Tel. 2510001/4. 


FUNDAÇÕES 


— Construções Técnicas, Lda. 
Pr. do Município, 13-3.º — Tel. 36 6506 


— Empresa de Sondagens e Fundações 
Teixeira Duarte, Ld.' 
Av. da República, 42-8.º 
— Telef. 77 70 86/7 e 77 12 65 — Lisboa 


— Empresa Ultramarina de Sondagens e 
Fundações 
R. S. Pedro de Alcantara, 1 — Tel. 324693 e 
213256 — Lisboa — Luanda — Lourenço Mar- 
ques — Beira, 


— Fundações Franki, Lda. 
R. Braamcamp, 13-3.º — Lisboa — Tel. 5341 12. 


— Johann Keller 
R. dos Bem Lembrados - Cascais —Tel, 28 o1 46 


— Sendagens e Fundações, A. Cavaco, Lda. 
R. Rodrigo da Fonseca, 62-r/c — Lisboa 1 
Tel. 5611 71/72/73. 
— Sondagens Ródio 
Rua de S. Bento, 644-3.º — Lisboa 
Tel. 6880 96/7/8 
— Sopecate 
R. do Arsenal, 146, 2.º — Lisboa. 
Tel. 340 10-320208 


IMPERMEABILIZAÇÃO 


— Empresa de Produtos Asfálticos 
R. Filipe Folque, 10-1.º — Lisboa — Tel. 7301 56. 


INSTRUMENTOS DE PRECISÃO 


— Nucleon 
Av. António Augusto de Aguiar, 165, r/c.-Dt. 
Tel. sg0r 05 — Lisboa 


— Wild Portugal, Ld.' 
Praça das Águas Livres, 8, s/| 6 — Lisboa 2 — 
Tel, 6891 12. 


ISOLAMENTO TÉRMICO E ACÚSTICO 


— Isola 
Av. António Augusto de Aguiar, 17-2., e 3.º Lx, 


— SETH, Ld.' 
R. Filipe Folque, 10-1.º — Lisboa — Tel. 7301 56 


SONDAGENS 


— Construções Técnicas, Lda. 
Pr. do Município, 13-3.º — Tel. 36 65 06 
— Empresa de Sondagens e Fundações Tei- 
xeira Duarte, Ld.' 
Avenida da República, 42-8.º 
Telef. 77 70 86/7 e 77 12 65 — Lisboa 


— Empresa Ultramarina de Sondagens e 


Fundações 
R. S. Pedro de Alcântara, 1 — Tel. 324693 e 
21 13 56 — Lisboa — Luanda — Lourenço Mar- 


ques — Beira, 


TECNICA— XXIX 


— Fundações Franki, Lda, 
R. Braamcamp, 13-3.º — Lisboa — Tel, 53 41 Ia 
a 


— Johann Keller 
R. dos Bem Lembrados - Cascais — Tel, 2801 46 


— Sondagens e Fundações, A. Cavaco, Lda. 
R. Rodrigo da Fonseca, 62-r/c — Lisboa 1 


Tel. 56 1171/72/73. 


— Sondagens Ródio, Ld.' 
Rua de S. Bento, 664-3.º — Lisboa 
Tel. 6880 96/7/8. 


— Sopecate 
R. do Arsenal, 146-2.º — Lisboa, 
Tela 40 10 — 320208 


CONSTRUÇÕES METALO-MECANICAS 


CONSTRUTORES METALO-MECÁÃ- 
NICOS 

— COMETNA 
R. Academia das Ciências, 5 
Tel. 3200 11 — Lisboa 

— Construções Metalo-mecânicas MAGUE 
Ld.'— Alverca 


— Sociedade Industrial Metalúrgica 
R. de S. Tiago, 1 


EQUIPAMENTO ELECTRICO 


— AEG — Lusitana de Electricidade 
R. João Saraiva, 4-Tel. 7r II 71/71 40 o1-Lisboa 


— ASEA 
Estrada do Prior Velho — Tel. 251 xr or 
— Sacavém 
R. do Campo Alegre, 144 —Porto — Tel. 62 106/9 


— Electrotécnicos Reunidos 
Avenida João XXI, 64 — Lisboa. 


— FRAPIL 
Cais de São Roque — Aveiro — Tel. 23071 


— J. F. de Azevedo e Silva & C.', Ld.' 
Rua S. Francisco Sales, 2, 1.º — Tel. 654165 
— Lisboa. 


— Hidrel 
R. do Alecrim, 47-A-B 
Tels. 325700 - 35966 — Lisboa 


— Jayme da Costa, Ld.' 
R. dos Correeiros, 14 — Lisboa — Tel. 327035 
Praça da Batalha, 12 — Porto — Tel, 2 28ra. 
Cx. Postal 5174 — Luanda — Tel. 4743. 


— Philips Portuguesa 
R. Joaquim António de Aguiar, 66 Lisboa 


TECNICA — XXX 


— Siemens — Companhia de Electricidade 
Avenida Almirante Reis, 65] 
Telef. 5369 21/0 — Lisboa 

— Soc. de Elec. Brown Boveri 
R. Sá da Bandeira, 481-2.º — Tel. 234 11— Porte 


— Sociedade Michãelis de Vasconcelos, Ld. 
— QOerlikon 


Av. Marquês de Tomar, 94 — Lisboa.? 
P. da Liberdade, 114 — Porto. 


— União Eléctrica Portuguesa 


R. Alexandre Herculano, 244 — Porto; 
R. Rosa Araújo, 35 — Lisboa, 


DISTRIBUIÇÃO DE ENERGIA 


— União Eléctrica Portuguesa 


R. Alexandre Herculano, 244 — Portos 
R. Rosa Araújo, 35 — Lisboa. 


ELECTRODOS 


— Electro-Arco, Ld. 


R. Silva Carvalho. 220 — Lisboa — Tel, 68 3649. 
R. do Bolhão, z16 — Porto — Tel. 2 12 77. 7a 


INSTALAÇÕES 


— ASEA 
Estrada do Prior Velho — Tel. 251 II OI 
— Sacavém 
R. do Campo Alegre, 144 —Porto —Tel. 62 106/9 


— |solux, Lda. 


R. Gonçalves Crespo, 33, 3.º — Lisboa 
R. de Sá da Bandeira, 551, 4.º — Porto, 


— JJ. F. de Azevedo e Siiva & C.”, Ld. 


Rva S. Francisco Sales, 2, '1.º — Tel, 654165 
— Lisboa. 


— J. Roma 
Praça da Figueira, 86'51/2/3 — Lx. 


— Jayme da Costa, Ld.' 
R. dos Correeiros, 14 — Lisboa — Tel. 327035. 
Praça da Batalha, 12 — Porto — Tel, 228 12, 
Cx. Postal 5174 — Luanda — Tel, 4743. 


— Siemens — Companhia de Electricidade 
Avenida Almirante Reis, 65 
Telef. 65 41 65 — Lisboa. 


— Standard Electrica 
Av. da Índia — Tel. 63 8171 — Lx. 


ISOLAMENTOS 


— Empresa Electro-Cerâmica, S.A.R.L. 
L. Barão de Quintela. 


SOLDADURAS 


— Electro Arco 


R. Silva Carvalho, 239 — Lisboa 
R. do Bolhão, 216 — Porto 


— J. Roma 
Praça da Figueira, 86/51/2/3 — Lx. 


TELECOMUNICAÇÕES 


— Empresa Técnics de Equipamentos Eléc- 
ricos, S.A.R.L. | 
R. Rodrigo da Fonseca, IIo, r/c — Lisboa — 
Tel. 6860 7a. 
k. Sá da Bandeira, 766-1.4E — Porto — Tel. 24818 


— Standard Electrica 
Av. Índia — Tel. 6381 71 — Lx. 


EQUIPAMENTO INDUSTRIAL 


— COMETNA 
Rua Academia das Ciências, 5 
Tel. 3200 11 — Lisboa 


— Frapil 
Cais de São Roque — Aveiro — Tel. 23071 


— DEMAG — Soc. Comercial Romar 
Trav. da Galé, q — Tel, 63 5670 — Lisboa. 
Santa Iria de Azoia 


— Guedes & Almeida, Lda, 
R. Aurea, 181-2.º — Tel. 3278 45 — Lisboa 


— Hidrel 
R. do Alecrim, 47-A-B 
Tels. 325700-3 59 66 — Lisboa 


— Kendal & Pinto Basto & C.”, Lda. 
R. Nova da Alfândega, 7-12 - Tel. 284 21 — Porto 
Av. 24 de Julho, 1 — Tel, 36 15 81 — Lisboa 


— Construções Metalo-Mecânicas Mague, 
Lda. 
Alverca do Ribatejo 


CORRENTES 


— Auto - Lusitania 
Av. da Liberdade, 73 -79 — Lisboa. 


— Harker Sumner 
L. do Corpo Santo, 18 
Tel. 32 4823 — Lisboa, 


— |JNA — Rolamentos de Agulhas 
Schaeffler, Lda — R. Sá da Bandeira, 740 — 
Porto. 
Tel: ar15 7/8/9. 


INDÚSTRIAS QUÍMICAS 


FELTROS 


— FANMAFEL —Fábrica Nacional” de Faltros 
Industriais 


Estrada de S. João — Ovar 


TEXTEIS 


— FANAFEL — Fábrica Nacional de Feltros 
Industriais, Ld.* 


Estrada de S. João — Ovar. 


LIVROS TÉCNICOS 


— TÉCNICA 
Avenida Rovisco Pais — Lisboa 
Tel, 77 54:49 


MATERIAIS DE CONSTRUÇÃO 


AÇOS ESPECIAIS 


— A. Johnson & €C.' (Portugal), Ld.* 
P. José Fontana, r1-1.º — Lisboa — Tel. 5370 15/ 
[6/7]8. 
R. Dr. António Granjo, 160/168-Porto-Tel. 5 46 68 
— Companhia Portuguesa de Trefilaria 


R. do Instituto Industrial, 18-1.º, Dt.” — Lisboa 
Tel. 67 12 24/5. 


— Heliaço 
R. S. Julião, 41-4.º — Tel, 36 78 17 — Lisboa. 


— PRISMA 
R. D* Luis 1, 20 Lisboa Tel. 6707 76 


BETÃO 


— Companhia Portuguesa de Trefilsria 
R. do Instituto Industrial, 18-1.º, Dt.* — Lisboa 


Tel. 67 12 24/5. 
Lisboa 1 


— Empresa Cimentos de Leiria 
R. Braamcamp, 7 — Lisboa. — Tel. s'gr 61/66. 


— Heliaço 

R.S. Julião, 41-4.º — Tel. 36 78 17 — Lisboa, 
— PRISMA 

R. D. Luis I, 20 Lisboa. Tel. 6701 76 
— Sociedade Portuguesa CAVAN 


R. de D. Estefânia, 94-A — Lisboa — Tel. 478 ta 


TECNICA — XXXI 


CIMENTOS 


— AGUIAR E MELO, LDA. 
Praç. do Município, 13-1.º — Tel. 821151/2 Lx 


— CIBRA — Cimentos Brancos 
Tv. do Corpo Santo, 15 — Lisboa — Tel. 3204 64. 
R. do Bonjardim, 205 — Porto — Tel. 25779. 


— Empresa de Cimentos de Leiria 
R. Braamcamp, 7 — Lisboa — Tel. 591 61/66. 


— Secil 
R. do Comércio, 156—Lisboa 2 —Tel. 328201/2/3 


MATERIAIS PRÉ ESFORÇADOS 


— CERVAL— Cerâmica de Vale de Lobos, Ld.* 
Av. Joaquim António de Aguiar, 41-5.º D. 
Tel. 53 11 38 — Lisboa. 


TECNICA — XXXII 


— Soc. de Pré-fabricação e Obras Gerais 
Novobra, Ld. 


Av. E. U. da América, 100-5.º-E. 
Tel. 77 48 32/ 76 64 46 — Lisboa. 


ARTIGOS DE ESCRITÓRIO E ESCOLAR 


MATERIAL DE DESENHO 


— J.J. Silva & Irmão, Lda. 
Rua Silva Carvalho, 232-B — Lisboa 
Tel. 652750 /6806 71 
— Nucleon 
Av. António Augusto de Aguiar, 165, r/c.-Dt.. 


Tel. 5301 05 — Lisboa 


— Wild Portugal 


Praça das Águas Livres, 8-s/1. 
Tel. 98 gr12 — Lisboa. 


Empresa de Sondagens e Fundações 
TEIXEIRA DUARTE, L.” 


AV. DA REPÚBLICA, 42, 8.º 
TELEFS. 777086/7 o = o L ! S B O PAN 


SONDAGENS GEOLOGICAS 
CAPTAÇÃO DE AGUAS SUBTERRANEAS 
CONSOLIDAÇAO DO SOLO 
INJECÇÕES DE CIMENTO 


FUNDAÇÕES DE TODOS OS TIPOS 
CAYES E TUNEIS 


BARRAGENS é PORTOS 


(UM QUARTO DE SÉCULO DE ESPECIALIZAÇÃO TÉCNICA) 


NON 
BOVERI 


Pára-raios 


Contadores de descargas 


Os vossos pára-raios funcionaram ? 


Quanias sobrelensões eliminaram ? 


O nosso contador de descargas 


dará resposta a estas perguntas 


Para cada problema de protecção, os nossos Serviços 


Técnicos podem indicar a solução mais apropriada 


soc. ve eLecr. BROWN BOVERI, coa. 


eee 7 lj) À DE SA DA BANDEIRA - 481-22 TEL. 2341414- PORTO 


